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I .  Introduction 

Adaptive  interference  estimation  and  suppression  algorithms 
have  recently  been  demonstrated  to  be  an  effective  method  of 
improving  the  performance  of  direct  sequence  spread  spectrum 
communication  systems  in  the  presence  of  narrowband  inter¬ 
ference  or  jamming.  Conventional  direct  sequence  spread 
spectrum  communication  receivers  rely  on  the  processing  gain 
inherent  in  this  communication  technique  to  provide  an  ad¬ 
vantage  over  the  intentional  jammer  or  other  co-channel  users. 
However,  Hsu  and  Giordano  [1],  and  Proakis  and  Ketchum  [2] 
have  demonstrated  that  by  taking  advantage  of  the  white  noise¬ 
like  properties  of  tne  direct  sequence  spread  spectrum  signal, 
any  of  a  number  of  spectral  estimation  algorithms  can  be 
employed  to  estimate  and  suppress  interference  which  is  narrow- 
band  relative  to  the  spread  spectrum  signal.  By  doing  this,  a 
significant  additional  advantage  over  the  interferer  is  gained. 

These  estimation  and  suppression  algorithms  have  in  common 
the  property  that  their  output  is  a  set  of  filter  coefficients 
for  a  transversal,  or  finite  impulse  response  (FIR),  filter 
which  processes  the  received  signal  prior  to  sequence  cor¬ 
relation,  as  shown  in  Figure  1.  The  estimation  and  suppression 
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r(t) 


Figure  1 

algorithms  attempt  to  design  a  transversal  filter  which  places 
a  notch  in  the  band(s)  occupied  by  the  narrowband  inter¬ 
ference.  Thus  these  filters  have  the  effect  of  significantly 
attenuating  the  interference  power,  as  well  as  introducing 
distortion  in  the  desired  signal. 

The  algorithms  considered  in  [2]  fell  into  two  classes: 
parametric  and  non-parametric.  The  non-parametric  algorithms 
make  use  of  classical  frequency  domain  techniques,  such  as 
described  in  [3]  to  provide  an  estimate  of  the  interference 
spectrum.  This  estimated  spectrum  is  then  used  to  compute 
the  coefficients  of  a  linear  phase  transversal  filter  which 
is  the  interference  suppression  filter. 

The  parametric  algorithms  considered  in  [1]  and  [2]  are 
based  on  modeling  the  interference  as  the  output  of  an  all¬ 
pole  filter  which  is  driven  by  a  white  noise  sequence.  The 
pole  positions  are  estimated  using  linear  prediction,  and  the 
resulting  pole  positions  are  used  as  the  zero  locations  of 
the  FIR  suppression  filter.  In  this  way,  the  intermediate 
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step  of  computing  a  spectral  estimate  is  avoided,  the  suppres¬ 
sion  filter  coefficients  being  computed  directly  from  the  data. 
In  view  of  their  greater  simplicity,  and  the  fact  that  they 
provide  generally  superior  performance,  the  parametric  methods 
only  will  be  considered  in  this  report. 

Both  of  these  methods  of  determining  a  suitable  suppression 
filter  were  considered  in  detail  in  [2],  and  will  not  be  con¬ 
sidered  further  here.  The  performance  of  these  suppression 
filters  was  given  in  terms  of  the  SNR  improvement  factor, 
which  is  the  ratio  of  the  SNR  with  the  suppression  filter  to 
the  SNR  without  the  filter.  No  explicit  error  rate  data  was 
given,  except  for  the  output  of  a  simulation  of  a  receiver 
which  uses  a  decision  feedback  equalizer  (DFE)  as  well  as 
suppression  filtering  in  a  multipath  environment.  Another 
question  which  was  not  considered  in  [2]  was  the  structure  of 
the  maximum  likelihood  receiver,  assuming  that  some  sort  of 
discrete  time  interference  suppression,  or  noise  whitening, 
filter  is  used.  It  is  these  two  subjects,  the  maximum  likeli¬ 
hood  receiver,  and  error  rate  performance,  which  are  the-  subject 
of  this  report.  In  addition,  some  corrections  to  results  con¬ 
tained  in  [2]  are  given  in  the  appendices  of  this  report.  The 
errors  involve  a  change  in  the  improvement  factor  in  the  case 
when  the  suppression  filter  is  linear  phase.  However,  they  in 
no  way  alter  the  conclusion  that  adaptive  interference  cancel¬ 
lation  is  an  effective  way  of  improving  the  performance  of 
direct  sequence  spread  spectrum  receivers  in  the  presence  of 
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narrowband  interference. 

Before  proceeding  with  the  body  of  this  report,  we  make 
some  definitions  which  are  useful  in  the  derivations  which 
follow. 

Some  Definitions 

Throughout  this  report,  equivalent  lowpass  representation 
of  bandpass  signals  and  systems  is  used.  Thus  the  bandpass 
signal,  y(t),  is  represented  by  the  lowpass  equivalent  signal, 
x(t),  satisfying  the  following  relation: 

j2irf  t 

y (t)  =  Re[x(t)  e  ] 

If  y (t)  is  a  deterministic  signal,  then  its  energy  is  given 
by: 


f 00  .00 

F-y  ■  J  y2(t)  dt  -  I  |  |*(t)|2  dt 
-  00  -  00 

If  y(t)  is  a  wide  sense  stationary  bandpass  process,  then  we 
can  define  its  autocorrelation  function  as: 

$yy(t)  =  E[y(t)  y (t  +  r)] 

The  autocorrelation  function  of  the  equivalent  lowpass  process 
is 


<J>xx(t)  "  \  x(t  +  T)  ] 


i 

i 


s 


These  two  autocorrelation  functions  are  related  by 

j27rfct 

<j>yy(t)  =  Re[<J>xX(t)  e  ] 

The  power  of  the  bandpass  process  is 

P  *  $  CO)  =  4>  (0) 

y  yy  xxv 

If  y(t)  is  a  sinusoid  in  the  band  of  interest,  of  amplitude 
A,  and  with  random  phase  <J> ,  then  y(t)  can  be  oppressed 

y(t)  =  A  cos[2ir(fc  +  f^)  t  +  4>] 

and  the  lowpass  equivalent  signal  is 

x(t)  =  A  exp{  j  [2tt  (f^)  t  +  <J> ] } 

The  respective  autocorrelation  functions  of  y(t)  and  x(t) 
are 


<J>yy(f)  =  V  cosI2lI(fc  + 


4>xx(t)  =  exp  [j  2  it  (f^)  x] 


and 
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The  power  is 


Assume  that  a  binary  baseband  communications  system  using 
antipodal  signals  with  pulse  waveform  p(t)  is  corrupted  by  a 
single  sinusoid  A  cos[2ir  f^t  +  <J>)  .  The  receiver  consists  of 
a  filter  matched  to  p(t).  The  signal  component  at  the  sampling 
instant  at  the  output  of  the  matched  filter  is: 


|  p (t  -  kxc)  p(t  -  ktc)  dr  =  E 


The  power  of  the  interference  component  is 


Tc  F'c  *ii(0)  =  T"  Tc  E 


Thus,  the.  signal -to- interference  ratio  is 
2F./A2  Tc 


Alternately,  assume  that  a  binary  bandpass  communications 

A 

system  uses  a  signaling  pulse  whose  lowpass  equivalent  is  p(t), 
and  is  corrupted  by  a  single  sinusoid  A  cos[2m(fc  +  f^)  t  +  <J>] 
whose  lowpass  equivalent  is  A  exp  [j  (2m  f^t  +  4>)}.  The  signal 
component  at  the  output  of  the  matched  filter  is 
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|  p*(x  -  kTc)  p(x  -  kxc)  dx  =  2E 
-  00 

and  the  power  of  the  interference  component  is 
2xc  E<j>ii(0)  =  A2  xc  E 
The  signal- to- interference  ratio  is 
4E/A2  x c 

Thus,  for  interference  with  the  same  power,  the  coherent 
bandpass  communications  system  achieves  twice  the  signal-to- 
interference  ratio  as  the  baseband  system. 

We  will  define  the  signal- to- interference  ratio  (SIR) 
of  a  bandpass  communications  system  as  the  ratio  of  signal 
power  to  equivalent  baseband  interference  power  in  a  baseband 
communications  system.  In  the  case  of  sinusoidal  interference, 
this  ratio  is 

SIR  =  4E/A2  xc 

1 1 .  Maximum  Likelihood  Receiver 

In  this  section,  we  consider  the  structure  of  the  maximum 
likelihood  receiver  for  detection  of  direct  sequence  spread 
spectrum  signals  in  wideband  (white)  Gaussian  noise,  and 
narrowband  interference.  In  its  most  general  form,  this 
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problem  has  no  straightforward  or  easily  implemented  solution. 
Here,  however,  we  make  the  simplifying  assumption  that  the 
receiver  contains  a  noise-whitening  filter  implemented  with  a 
tapped  delay  line,  and  whose  coefficients  are  computed  using 
one  of  the  parametric  methods  given  in  [2J.  The  receiver 
structure  for  transmission  over  a  nondispersive  channel  will 
be  considered  first,  then  these  results  will  be  generalized 
to  include  channels  with  multipath  distortion. 


Nondisversive  Channel 


The  model  used  for  the  transmitter  and  nondispersive 
channel  is  shown  in  Figure  2.  The  input  to  the  transmitter 
is  a  binary  source  which  generates  symbols  u^ ,  chosen  with 
equal  probability  from  the  binary  alphabet  (-1,  +1).  The 
source  sequence  Uj  is  multiplied  by  the  modulating  PN 
sequence,  c^,  and  the  product  modulates  an  impulse  generator. 
The  subscript  j  of  the  source  sequence  is  related  to  the  sub¬ 
script  k  of  the  PN  sequence  by:  j  =  (k/Lj ,  where  L  is  the 


Figure  2 
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processing  gain  expressed  in  chips  per  bit,  and  JxJ  denotes 
the  greatest  integer  not  greater  than  x. 

The  output  of  the  impulse  generator  is  shaped  by  the 
filter  P(f),  and  is  transmitted  over  the  channel  which  is 
assumed  to  have  a  flat  frequency  response  and  attenuation  a. 

i 

The  transmit  filter  P(f)  is  assumed  to  be  such  that  it  intro¬ 
duces  no  intersymbol  interference  between  the  elementary  chip 
pulses  transmitted  on  the  channel  (interchip  interference) . 

That  is,  if  P(f)  is  followed  immediately  by  its  matched 
filter  P*(f),  whose  impulse  response  is  p*(-t),  and  the  output 
of  the  matched  filter  is  sampled  at  the  appropriate  moment, 
then  energy  from  only  one  chip  pulse  effects  the  sample  value. 

This  condition  can  be  written  as: 

p(x)  p*(x  -  kxc)  dx  =  2Ec  5ko  (1) 

-  00 

where  x  is  the  reciprocal  of  the  chip  rate,  E  is  the  energy 
c  ^ 

in  the  bandpass  chip  pulse,  and  6ko  is  the  Kronicker  delta: 

6,  =  I5  2  •  Finally,  we  assume  that  the  receiver  has 
ko  10,  k  f  0 

perfect  knowledge  of  the  phase  of  the  received  signal,  so  that 
for  the  purposes  of  this  discussion,  there  is  zero  phase  shift 
in  the  channel,  and  a  is  a  real  constant. 

White  Gaussian  noise  and  narrowband  interference  are 
added  to  the  signal  at  the  receiver.  The  double-sided  power 
spectral  density  of  the  bandpass  white  noise  process  is  NQ/2, 
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and  the  lowpass  equivalent  process,  n(t),  has  double-sided 
power  spectral  density  NQ,  with  corresponding  autocorrelation 
function: 


4»nn  (t  )  =  E [n (t)  n*  (t  -  t)]  =  Nq  6(x)  (2) 

The  lowpass  equivalent  narrowband  interference  has  autocor¬ 
relation  function: 


♦  u(t)  «  \  E[i(t)  i*(t  -  x)] 


(3) 


The  received  signal  during  the  symbol  epoch  is  given 

by: 


Cj+llL-l 

r^(t)  =  a  I  U(k/Lj  ck  p(t  “  kTc)  +  n(t)  +  i(t) 

JL  (4) 


For  the  purposes  of  deriving  the  maximum  likelihood 
receiver  for  the  transmitted  message  we  will  assume  that  prior 
to  further  processing,  the  receiver  filters  the  received 
signal  with  a  tapped  delay  line  whitening  filter  whose  impulse 
response  is: 

P 

h(x)  »  l  h (m)  5 (x  -  mx  )  (5) 

m=0  c 

It  will  be  demonstrated  subsequently  that  this  is  equivalent 
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to  a  receiver  which  contains  a  sampler  with  sampling  period 
t c ,  followed  by  a  sampled  data  FIR  filter  with  coefficients 
hm>  m  =  0 »  • • • ,  P. 

The  output  of  the  whitening  filter  during  the  j symbol 
epoch  is  given  by: 

P 

W.(t)  =  l  h (m) 

J  m=0 

Cj+1)L-1  P 

=  a  l  l  h  (m)  u,v/T|  Ci  p[t  -  (m  +  k)  x  1 

k+jL  m=0  K  c 


P 

+  l  [n(t  -  mt  )  +  i  (t  -  mx  )]  h(m)  (6) 
m=0  c  c 

It  is  this  signal,  Wj(t),  for  all  j  over  the  message  duration, 
to  which  the  maximum  likelihood  criterion  will  be  applied  in 
order  to  determine  the  structure  of  the  maximum  likelihood 
receiver. 

Before  proceeding,  we  will  make  the  following  definitions 

x(t)  =  p(t)*h(t)  (7) 


6(x  - 


mxc) 


r(t  -  x)  dx 


and 


z(t)  =  [n(t)  +  i(t)]*h(t) 


(8) 
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where  the  operation  *  indicates  convolution.  The  signals 
x(t)  and  z(t)  are  the  components  at  the  output  of  the  filter 
defined  in  (5)  due  to  the  transmitted  signal,  and  the  noise 
and  interference,  respectively.  Furthermore,  without  loss 
of  generality,  we  will  assume  that  a  given  message  has  dura¬ 
tion  -NT  <  t  <  NT,  or  -NLtc  <  t  <_  NLtc  -  1,  where  T  =  Lx 
is  the  symbol,  or  bit,  duration.  Thus,  the  received  signal 
is  written: 

NL-1 

r(t)  =  a£  u  ...  c.  p(t  -  kx  )  +  n(t)  +  i  (t)  (9) 

k=-NL  K/L  K  c 


and  at  the  output  of  the  whitening  filter,  the  signal  is: 


NL-1 

w(t)  =  al  u  k/,L  x ( t  -  kxc)  +  z(t)  (10) 

k=  -  NL 

In  order  to  obtain  the  likelihood  function,  p[w(t)|uj, 
which  must  be  maximized  in  order  to  minimize  the  probability 
of  error,  w(t)  will  first  be  expanded  using  the  Karhunen-Lobve 
series  expansion  as  described,  for  example,  in  [3].  This 
process  will  provide  a  set  of  signal  samples,  w^  which  are 
uncorrelated.  We  assume  that  z(t)  is  a  Gaussian  process;  thus, 
the  w^  are  independent  Gaussian  random  variables.  The  Karhunen- 
Lofeve  expansion  for  w(t)  is  given  by: 

N 

w(t)  «  lim  l  w,  f  (t) 

N-*°°  k=  1 


(ID 
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where 


NL-1 

“k '  e“L  ui*/y  c«  xk,i  *  zk 


(12) 


The  functions  £ ^ ( t: )  are  an  orthonormal  set  which  satisfy 


(NL-l)t, 


-NLt, 


fk(T)  *zz(t  •  T)  dT  =  h  fk(t) 


and  z ^  and  x ^  ^  are  given  by: 


(NL- 1) t 


•NLt, 


z(t)  fR*(t)  dt 


and 


lk,£ 


(NL-1)  r 


C  x(t  -  £tc)  fk*(t)  dt 


-NLt, 


i. 

Thus  the  likelihood  function  for  the  p  possible  transmitted 
sequence.  u^,  can  be  written: 


N 


P Cw !  U  )  =  lim  [  n  2  it  A,] 

V  _  i  K 


-1/2 


N+co  k=l 


,  N  NL-1  2 

• exp[-  Ik!,vk-  t.^LUP.  t/L  c*  xk.«l  /Xk 


(13) 
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If  we  assume  that  the  whitening  filter,  h(t),  has  done  its  job 
perfectly,  so  that  z(t)  is  a  white  Gaussian  process,  it  can  be 
shown  that 


p(w|up) 


NL-1 

uak u  m  c*  x(t 


dt] 


(14) 


Maximizing  (14)  over  all  possible  transmitted  sequences,  u^, 
can  be  accomplished  by  maximizing  £n[p(w|Up)]  which  is  equi¬ 
valent  to  minimizing  the  metric 


NL-1 

=  2Re [  l 

k=-NL 


u 


*  tk/Lj 


w(t)  x*(t 


kxc) 


dt  ] 


NL-1  NL-1 

k=-NL  £=-NL  UP»Lk/LJ  P»  ll/l\ 


c 


Z 


x*(t  -  kTc)  x(t  -  Jt  t  c )  dt 


(15) 


Note  that  since  both  Uj  and  c^  are  chosen  from  a  binary  real 
alphabet,  (-1,  1),  the  complex  conjugate  indications  for  the 
u’s  and  c's  in  (15)  can  be  dropped. 

Next,  we  define  the  sample  value  y<K: 
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yk  =  Re  j  w(x)  x*(t  -  kxc)  dx 


NL-1 


-  «Z 

S,  =  -NL 


«.TC)  x*(x  -  kxc)  dx 


+ 


z(x)  x*(x  -  kxc)  dx 


’  Re[2oEc  “(c^n-mJ/Lj  ck*n-m  h(m)  h‘(n)l 


P  P  r» 

+  Re{  l  l  hm  hn*  [n(x  -  mxc)  +  i(x  -  mxc) ] 
m=0  n=0  ' _  v 


•  p*[x  -  (n  +  k)  x c ]  dx) 


P 

-  Re  [2otE  £ 
m=~  p 


uL(k-m)/Lj 


ck-m  s(m) 


+ 


[n(x  -  mxc)  +  i(x  -  mxc)] 


•  p* ( x  -  kxc)  dx] 


(16) 


where 
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g(m)  &  l  h (m  +  n)  h*(n),  m  =  -P,  P  (17) 

n=0 


The  sample  value  y^  can  be  obtained  by  passing  r(t)  through 
the  whitening  filter  h(x)  given  in  (5)  (see  Figure  3),  then 
through  its  matched  filter,  h*(-x),  and  through  the  chip  pulse 
matched  filter,  p*(-x).  The  output  of  this  filter  would  be 
sampled  at  rc  second  intervals,  and  the  real  part  of  the  sample 
values  taken,  thus  forming  y^.  However,  since  the  desired 
signal,  u|jc/IJ  ck*  *s  real»  there  is  an  alternate  strategy 
which  is  preferable,  shown  in  Figure  4.  The  first  filtering 
operation  to  be  performed  in  the  receiver  is  the  matched 
filter  operation,  /°°  r(r)  p*(x  -  kxc)  dx .  The  output  of  the 

-  oo 

matched  filter  is  sampled,  and  the  real  part  of  the  sampler 
output  is  taken.  Thus,  we  have  eliminated  the  need  to  perform 
spectral  estimation  on  a  complex  valued  signal  and  can  now  use 
a  discrete  time  whitening  filter  with  real  valued  tap  weights. 

Next,  noting  that 


P  P 

x*(t  -  kx  )  x ( t  -  Jtx  )  dt  =  l  l  h  (m)  h(n) 
>  c  m=0  n=0 


J  p*[t  -  (k  +  n)  xc)  pft  -  («,  +  m)  x c ]  dt 


2Ec  l  6£k+m 

m=  -  p 


(18) 
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we  can  write  the  likelihood  function  as: 


NL-1 


ap  ■ 


JNL  t2UP.  Lk/Lj  Ck  yk  '  2“Ec  J  Up,Lk/Lj 


P  (j*l)L-l 

'  4“Ec  Jl  JjL  Vi  UP.LCk-"»Aj  Ck  ck-n 


The  first  term  inside  the  brackets  in  (19)  indicates  that 
the  output  of  the  composite  whitening  filter-matched  filter  is 
processed  by  a  sequence  correlator  which  multiplies  by  the 
PN  sequence  c^,  and  sums  over  the  bit  interval  T.  If  there 
were  no  interference,  and  therefore  the  whitening  filter  was 
such  that  g(m)  =  0,  m  f  0,  then  this  operation  alone  would 
constitute  the  maximum  likelihood  receiver  and  the  output  of 
the  sequence  correlator  would  form  the  decision  variable  on 
which  bit  decisions  would  be  made,  on  a  bit-by-bit  basis. 
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However,  because  of  the  dispersion  caused  by  the  noise¬ 
whitening  filter,  depends  not  only  on  the  bit  currently 
being  transmitted,  but  also  on  adjacent  bits.  For  this 
reason,  it  is  necessary  to  consider  the  received  sequence  in 
its  entirety,  rather  than  on  a  bit-by-bit  basis,  when  deciding 
which  transmitted  sequence  was  most  likely  to  have  been  trans¬ 
mitted. 

The  Viterbi  algorithm  is  an  efficient  method  of  accom¬ 
plishing  this  maximum  likelihood  sequence  estimation  without 
waiting  for  the  entire  sequence  u  to  be  transmitted.  If  we 

wait  for  the  entire  sequence  of  2N  bits  to  be  transmitted,  and 

2N 

consider  the  metric,  A  ,  for  each  of  2  possible  transmitted 

r 

sequences,  the  result  is  a  very  complex  receiver  with  large 

decoding  delay.  The  Viterbi  algorithm  reduces  the  delay  and 

complexity  at  the  receiver  by  constructing  a  state  trellis 

through  which  the  receiver  moves  as  it  performs  the  maximum 

likelihood  detection.  The  trellis  has  2N  stages  with  2^P/^ 

nodes  at  each  stage,  where  P  is  the  memory,  in  chip  intervals, 

introduced  by  the  whitening  filter;  L  is  the  processing  gain; 

and  fx"!  denotes  the  least  integer  not  less  than  x.  At  each 

stage  in  the  trellis,  which  corresponds  to  processing  of  a 

single  received  bit,  the  receiver  computes  two  branch  metrics 

for  each  node.  Thus  at  each  stage,  the  receiver  considers 

2^/1!  +1  sequences,  and  (2N)2^>^+*  for  the  entire  sequence. 

t*H 

The  branch  metric  which  must  be  computed  for  the  P  branch 
at  the  jth 


stage  is: 
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X-‘  =  1-  2up,j  ck  ^k  -  2l':c”(up,j  Ck>‘  8(°) 


PJ 


( j  + Y  L- 1 
k= j  L 


(j+l)L-l  P 

-  4a!:,.  I  *  )  u  .  u  i,,  ,  I  c,  c,  g(m) 

C  k-jL  m=l  P*J  P,|(k-*)/y  k  k-m 

(20) 

A  block  diagram  of  the  receiver  which  results  from  the  maximum 
likelihood  criterion  is  shown  in  Figure  5,  as  well  as  a  state 
trellis  for  P/L  =  1  in  Figure  6. 

It  is  interesting  to  note  that  the  maximum  likelihood 
criterion  leads  to  a  branch  metric  which  does  nothing  to 
counter  the  effects  of  pulse  dispersion  on  the  chip  level. 

To  see  this,  assume  that  L  >>  P,  so  that  the  contribution  to 
the  rightmost  term  in  (20)  due  to  summands  where  |_(k  -  m)/Lj  f  j 
is  negligible.  If  this  is  the  case,  for  a  given  j  this  term  is 
constant  over  all  sequences,  and  for  this  reason  it  has  no 
effect  on  the  sequence  decision.  In  other  words,  the  inter¬ 

symbol  interference  is  negligible.  In  this  case,  the  branch 
metric  reduces  to: 


A  . 
PJ 


+ 


(j  +  lJL-1 

l  i  c,  yk 

k=j L  k  K 


(j+lH-l 
k= j  L 


2up,j  ckt2oEc  uj  ck  SWI 


u*m-i 

k=  j  L 


2u 


P.j  Ck[2QHc  mJ_p  uL(k-m)/Lj  Ck-m  g(m)] 
m^O 


+  noise  term 


(21) 
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The  first  term  in  (21)  is  the  desired  signal,  while  the  second 
term  is  a  self-noise  or  "interchip"  interference  term  which  is 
always  present  whether  or  not  there  is  intersymbol  interference. 
This  self-noise  term  can  seriously  degrade,  and  in  many  cases 
turns  out  to  limit,  the  performance  of  maximum  likelihood 
detection. 

In  a  sense,  the  fact  that  the  ML  criterion  does  nothing 
to  mitigate  "interchip"  interference  is  due  to  the  fact  that 
over  a  given  bit  interval,  there  is  only  one  possible  chip 
sequence  (assuming  that  the  receiver  is  synchronized  to  the 
transmitted  PN  sequence).  Once  the  receiver  has  considered 
the  two  possible  symbol  values  for  that  interval,  it  has  ex¬ 
tracted  all  available  information  about  the  channel  from  the 
received  sequence.  This  view  leads  to  the  speculation  that 
it  may  be  desirable  to  introduce  coding  that  will  increase 
the  number  of  possible  sequences  during  a  given  bit  interval. 

Another  observation  which  can  be  made  about  (21)  is  that 
this  metric  is  the  same  metric  which  is  obtained  by  applying 
the  maximum  likelihood  criterion  to  transmission  over  a  dis¬ 
persive  channel  corrupted  by  white  noise  at  the  receiver  input 
when  the  channel  is  modeled  as  a  tapped  delay  line  with  weights 
h(m).  In  other  words,  for  the  purposes  of  digital  transmission 
and  maximum  likelihood  detection,  a  dispersive  channel  cor¬ 
rupted  by  white  noise  is  equivalent  to  a  nondispersive  channel 
corrupted  by  colored  noise  whose  power  spectral  density  is  the 
inverse  of  the  magnitude  squared  frequency  response  of  the 
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equivalent  dispersive  channel. 

Dispersive  Channel 

To  obtain  the  maximum  likelihood  receiver  for  the  dis¬ 
persive  channel,  corrupted  by  white  noise  and  narrowband 
interference,  we  replace  p(t)  by  p'(t)*c(t),  where  p'(t)  is 
the  transmitted  pulse,  and  c(t)  is  the  channel  impulse  response. 
The  impulse  response  c(t)  is  ip  general  a  complex  valued 
function,  and  is  also  possibly  time- varying .  In  this  case 
x(t)  and  r(t)  are  given  by: 

x(t)  *  p ’ (t) *c(t) *h(t)  (22) 


and 


r(t) 


NL-1  r°° 

IcJnl  UlyLJ  Ck  J  P'(T)  ctt  •  T  •  kTc) 

-  oo 

(23) 


The  input  to  the  sequence  correlator,  y^,  is  given  by: 

L 

yk  =  2Ec  mJ_L  uJ_(k-m)/Lj  ck-m  Rxx(m) 


P 

+  Re  I 
m=-p 


[n(t 


mxc)  +  i(t  -  mtc)] 


•  p* (t  -  ktc)  dx 


(24) 
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where 


Rxx00  -  j  X(T)  x*(T  -  kxc)  dr  (25) 

Note  that  since  c(t)  is  in  general  complex,  we  can  no  longer 
take  the  real  part  of  the  output  of  the  chip  pulse  matched 
filter.  If  we  follow  the  chip  pulse  matched  filter  with  the 
channel  matched  filter,  whose  impulse  response  is  c*(-nTc), 
then  we  can  take  the  real  part  of  the  channel  matched  filter 
output,  and  can  still  use  an  interference  suppression  filter 
with  real  valued  tap  weights.  This  receiver  structure  is 
shown  in  Figure  7. 

Placing  the  channel  matched  filter  prior  to  the  noise¬ 
whitening  filter  has  some  disadvantages,  however.  The  first 
is  that  since  we  no  longer  have  a  nondispersive  channel,  the 
component  of  the  received  signal  due  to  the  transmitted  signal 
no  longer  has  a  flat  power  spectral  density.  The  spectral 
estimation  algorithm  which  determines  the  noise-whitening 
filter  tap  weights  will  now  react  to  the  channel- induced 
distortion  of  the  signal  as  well  as  to  the  interference.  The 
channel- matched  filter,  having  the  same  magnitude  response  as 
the  channel,  will  reinforce  any  distortion  caused  by  the 
channel  and  thus  increase  the  undesired  response  of  the 
whitening  filter  to  the  channel.  The  second  disadvantage  of 
this  structure  is  that  if  the  channel  is  not  known  a  priori 
at  the  receiver  then  it  must  be  estimated.  A  straightforward 


Figure  8 
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way  of  doing  this  is  to  use  a  RAKE  (Figure  8)  matched  filter 
[4];  however,  it  is  not  practical  to  implement  the  RAKE  prior 
to  the  noise»whitening  filter  for  two  reasons.  The  first  is 
that  the  PN  sequence  correlation  is  performed  internally  to 
the  RAKE  so  that  its  output  is  the  narrowband  data  signal, 
and  it  is  no  longer  possible  to  perform  noise-whitening  on 
this  signal.  Even  if  it  were  possible  to  remove  the  sequence 
correlator  from  the  RAKE,  this  structure  would  still  be  un¬ 
desirable  because  the  large  interference  component  at  the 
RAKE  input  makes  the  job  of  channel  estimation  unnecessarily 
difficult. 

A  more  desirable  arrangement  is  shown  in  Figure  9.  Here 
the  chip  pulse  matched  filter  is  followed  by  a  sampler  and 
then  the  noise-whitening  filter.  The  noise-whitening  filter 
must  now  have  complex  valued  tap  weights  since  we  are  dealing 
with  complex  valued  signal  samples.  The  noise-whitening 
filter  is  followed  by  a  RAKE  matched  filter  which  adapts  itself 
to  have  the  response  of  a  filter  matched  to  the  cascade  of  the 
channel  and  the  noise-whitening  filter.  Although  it  may  seem 
unnecessary  to  have  the  RAKE  adapt  itself  to  the  whitening 
matched  filter  as  well  as  the  channel  matched  filter,  there  is 
no  other  straightforward  way  of  obtaining  filters  matched  to 
both  the  channel  and  the  whitening  filter.  The  RAKE  operates 
in  such  a  way  that  it  adapts  its  filter  response  to  match  the 
total  response  of  the  channel  and  receiver  between  the  trans¬ 
mitter  and  the  RAKE.  For  example,  if  the  whitening  filter 


were  followed  by  its  matched  filter,  and  then  by  the  RAKE, 
the  RAKE  would  adapt  itself  to  the  response  of  a  filter  matched 
to  the  channel,  the  whitening  filter,  and  its  matched  filter. 
For  this  reason,  we  must  allow  the  RAKE  to  match  the  channel 
and  the  whitening  filter.  Simulation  results,  which  will  be 
presented  in  a  subsequent  section  demonstrate  that  the  penalty 
for  implementing  the  noise-whitening  matched  filter  adaptively 
within  the  RAKE  in  this  way  is  negligible. 

III.  Bit  Error  Probability  for  Symbol-by- Symbol  Detection 

The  first  results  given  in  this  section  are  bit  error 
probability  expressions  for  symbol-by- symbol  detection  based 
on  a  Gaussian  assumption  about  the  total  noise,  including 
self-noise,  at  the  sequence  correlator  output.  Two  receivers 
are  considered:  the  conventional  direct  sequence  spread 
spectrum  receiver  shown  in  Figure  10  which  consists  of  a 
filter  matched  to  the  chip  pulse  and  a  sequence  correlator; 
and  an  enhanced  spread  spectrum  receiver  shown  in  Figure  11 
which  in  addition  to  the  chip  pulse  matched  filter  and  sequence 


correlator  also  has  an  interference  suppression  filter.  The 
interference  suppression  filter  consists  of  a  noise-whitening 
filter,  or  a  noise-whitening  filter  in  cascade  with  its 
matched  filter. 
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Receiver  with  No  Suppression  Filter 

The  decision  variable  at  the  output  of  the  sequence 
correlator  in  Figure  10  is  given  by: 


z  j  =  Re 


U  +  lJL-1 


f 


L  *=jL 


I  r. (t)  p* (t  -  £xc)  dx 


=  au . 
J 


Re 


‘(j  +  l)L-l 
.  £  =  j  L 


(j+r)L-i 
k= j  L 


p(x  -  kxc) 


•  p*(x  -  £xc)  dx 


Cj  +  1JL-1 

+  I  c,  [n(x)  +  i  (x)  ]  p*(x  -  kx  )  dx 

k-jl  L  . 


2aEc  LUj  +  Re 


(j  +  l^L-1 


k= j  L 


ck  |  [n (t)  +  i  (t)  ] 


p*(x  -  kx  )  dx 


(26) 


If  we  make  the  reasonable  assumption  that  z^  is  a  Gaussian 
random  variable,  then  the  mean  and  variance  of  (26)  are  suf¬ 
ficient  to  give  the  error  rate.  These  are: 


E(Zj)  =  2aEc  LUj 


(27a) 


i 

i 


B 
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var(z  • )  =  N.  +  I . 
3  3  J 


(27b) 


where  N.  =  2L  E  N  is  the  variance  of  the  white  noise  term 

J  CO 

at  the  output  of  the  sequence,  and 


I . 
3 


(j  +  r)L-l 

Jt= j  L 


wy1-1 

k= j  L 


Ck  C* 


.  oo 

2<p.  .  (t 

JJ  ii 


x) 


•  p*(t  -  krc)  p (t  -  £tc)  dt  dx  (28) 

is  the  variance  of  the  interference  term.  The  double  integral 
in  (28)  can  be  evaluated  first  by  noting  that  it  represents 
the  autocorrelation  of  the  interference  component  of  the  output 
of  the  chip  pulse  matched  filter,  p*(-x).  Denoting  this  signal 
component  by  i(t),  this  autocorrelation  function  can  be  written 

4>^(T)  -  f  |  P(f)  |  2  Si(f)  ej2nfz  df  (29) 

-  OO 

2 

where  |P(f)|  is  the  magnitude  squared  frequency  response  of 
the  chip  pulse  matched  filter,  and  S^(f)  is  the  power  spectral 
density  of  the  interference  at  the  input  of  the  chip  pulse 
matched  filter.  Assuming  that  the  chip  pulse  satisfies  the 
Nyquist  criterion  for  lack  of  intersymbol  interference,  and 
that  |  P(f)  j  2  -  constant ,  -  <  f  <  j~—  ,  then  |  P(f)  |  2  = 

2tc  Ec’  s0 
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foo 

Ct)  -  2t  E„  S. (f)  eJ 
11V  '  C  C  J  1 v  1 


2- irfx 


df 


=  2t  E„  <}> .  .  ( t ) 
c  c  yn v 


(30) 


Thus,  1^  is  given  by: 


I  .  =  2t  E„  L<|>.  .  (0) 
j  c  c  11^' 


(31) 


The  bit  error  probability  is  then  given  by: 


P u  =  Pr(z, 


(z.  <  0|u.  =  +1)  *  i  erfc  /y 
J  J  L 


(32) 


where 


[E(z.)]2  L  Ec  a2 

Y  "  2var (Zj )  "  N0  +  Tc  ^ii^ 


(33) 


Receiver  with  Suppression  Filter 

The  decision  variable  at  the  output  of  the  sequence 
correlator  in  Figure  11  when  the  excission  filter  consists  of 
the  noise-whitening  filter  plus  its  matched  filter  is: 
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•  P*[T  -  (k  -  £)  Tc]  dxj  (34) 

When  the  noise-whitening  filter  alone  is  used,  the  decision 
variable  is: 

cj*yL-i  p 

Zj  •  2“Ec  Lh(°)  uj  *  2“Ec  k=iL  u|_(k-«)/Lj 

•  ck  ck- £  h(£) 

(j  +  l)L-l  P  r°° 

+  l  l  c.  h(i)  Re  [n(x)  +  i  ( t )  ] 

k= j  L  £=0  K  J 

•  p*[t  -  (k  -  i)  x c ]  dx  (35) 

In  finding  the  mean  and  variance  of  z ^ ,  we  assume  that  the 
starting  phase  of  the  PN  sequence  is  a  uniformly  distributed 
random  variable,  which  is  uncorrelated  with  the  phase  of  the 
information  sequence,  so  that  the  first  PN  chip  in  any  given 
information  bit  occurs  with  equal  probability  at  any  position 
in  the  PN  sequence.  We  also  assume  that  the  autocorrelation 
of  the  PN  sequence  is  E[ck  c4]  =  «kr  Another  way  of  stating 
these  assumptions  is  that  the  channel  symbol  for  the  j1"*1  trans¬ 
mitted  bit  is  given  by  Uj  c y  where  u^  is  the  information  bit 
and  Cj  is  a  random  binary  vector  of  length  L  with  elements  chosen 
from  the  binary  alphabet  (-1,  1).  The  vector  c ^  can  be  any 
of  the  2^  possible  binary  vectors  of  length  L,  with  equal 
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probability . 

Thus,  in  the  whitening  filter  only  receiver,  the  mean 
and  variance  of  z^  are: 

E  (z  . )  =  2ctE  L  h(0)  u.  (36) 

J  C  J 

P  7 

var(z.)  =  2L  N  E_  f  [h  ( z 

y  0  c  *=o 

p  p 

+  2L  x  E  l  l  h(£)  h  (m)  <|>.  .[(m-A)  t  ] 

£=0  m=0 

p 

+  4a2  L  E2  l  [h  (£) ] 2  (37) 

£=1 

The  signal- to-noise  ratio,  defined  as  E  (z . ) /2var (z . )  is  then: 

j  j 

L  Ec  h2 (0)  a2 

Yj  p  p  p  *  p 

N  l  h2(£)+x  l  l  h(£)h(m)4>..[(m-£)x]+2a2E  l  h2(£) 
°  £=0  C  £=0  m=0  ii  C  £-1 

(38) 

and  the  bit  error  probability  is  again  given  by: 

Pb  =  \  erfc  /YT 

When  the  receiver  contains  the  filter  matched  to  the  whitening 
filter  as  well  as  the  whitening  filter,  the  mean  and  variance 
of  z.  are: 
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E(z.)  =  2otE  Lg (0)  u. 
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(39a) 


var (z . )  =  2L  N  E  £  gZ(£) 

J  £  =  -  p 


P  P 
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The  signal- to-noise  ratio  in  this  case  is 


y  i  u 


L  Ec  g2(0)  a2 
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(40) 


For  a  detailed  derivation  of  these  formulas,  see  Appendix  A. 

We  next  give  some  performance  results  obtained  by  use  of 
these  bit  error  probability  expressions.  In  all  of  these 
results,  we  have  taken  tc  =  1,  and  a  =  1.  The  results  are 
plotted  against  E^/NQ  to  aid  in  comparison  with  the  performance 
of  ideal  antipodal  signaling  without  interference.  The  inter¬ 
ference  is  a  sum  of  100  sinusoids  occupying  201  of  the  band, 
between  0  and  0.1  Hz.  The  results  are  given  for  processing 
gains  of  10,  20,  30,  and  60  chips  per  bit.  The  performance  of 
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the  receiver  without  interference  suppression  as  given  by 
(32)  and  (33)  is  not  plotted  since  at  an  SIR  of  -20dB  and 
the  processing  gains  given  above  this  probability  of  error 
is  not  much  below  0.5  for  any  E^/Ng. 

Figure  12  shows  the  performance  with  a  suppression  filter 
which  uses  a  4-tap  predictor,  and  no  matched  filter,  at  an  SIR 
of  -20dB/chip.  At  the  low  processing  gain  of  10  the  sup¬ 
pression  filter  has  provided  only  slight  improvement  over  the 
receiver  with  no  suppression;  at  a  processing  gain  of  60  there 
is  a  significant  improvement  over  the  receiver  without  sup¬ 
pression,  but  clearly  there  is  room  for  improvement  relative 
to  the  performance  without  interference. 

Next,  Figure  13  gives  the  performance  arrived  at  by  the 
Gaussian  assumption  for  the  same  interference  and  receiver 
as  that  in  Figure  12,  with  the  exception  that  now  the  receiver 
includes  a  filter  matched  to  the  suppression  filter.  This  has 
provided  a  marked  improvement  in  performance  over  that  obtained 
without  the  matched  filter,  particularly  at  the  higher  pro¬ 
cessing  gains.  It  is  notable,  however,  that  for  all  these 
processing  gains  the  bit  error  rate  always  bottoms  out  at  some 
finite  probability  of  error.  The  source  of  this  behavior  will 
be  addressed  presently. 

Figure  14  shows  the  performance  of  the  same  receiver  with 
lower  interference  power  resulting  in  an  SIR  of  -lOdB/chip. 

This  drop  in  interference  power  results  in  an  expected  per¬ 
formance  improvement.  The  performance  bottoms  out  at  lower 
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□  PROCESSING  GAIN*  iO 
O  PROCESSING  GAIN*  20 
*  .PROCESSING  GAIN*  30 
+  PROCESSING  GAIN-*  60 


Figure  12 


Figure  13 


Figure  14 


error  probabilities  at  all  processing  gains;  at  a  processing 
gain  of  60,  the  bottoming  out  is  off  the  scale  of  the  plot. 

_  7 

At  the  higher  processing  gain,  the  Ej3/Nq  requirement  for  10 
bit  error  probability  is  degraded  by  about  3dB  from  ideal 
binary  signalling,  compared  with  about  6dB  for  an  SIR  of 
- 20dB. 

The  4- tap  predictor  used  so  far  is  close  to  the  minimum 
order  filter  required  to  place  a  notch  in  the  band  occupied 
by  the  interference  and  thus  approximate  a  noise-whitening 
filter,  as  established  in  a  previous  report  [2].  It  was  also 
shown  that  by  increasing  the  filter  order  modest  gains  in  the 
SNR  improvement  factor  could  be  obtained  over  the  4- tap  pre¬ 
dictor.  In  the  frequency  domain,  the  effect  of  increasing 
the  predictor  order  is  a  notch  which  more  closely  approximates 
the  spectral  shape  and  location  of  the  interference.  Figure 
15  shows  the  error  probability  for  the  same  interference  as 
Figure  13,  except  now  the  predictor  is  a  15-tap  predictor. 

This  has  led  to  a  significant  gain  in  performance  over  Figure 
13;  in  fact,  the  performance  of  a  15-tap  predictor  in  -20dB 
of  interference  is  slightly  superior  to  that  of  a  4-tap  pre¬ 
dictor  in  -lOdB  of  interference. 

We  have  only  considered  here  the  performance  of  narrow- 
band  interference  suppression  when  the  interference  is  con¬ 
tained  within  a  single  band.  However,  it  was  demonstrated 
in  [2]  that  similar  values  of  the  SNR  improvement  factor, 
which  is  proportional  to  the  SNR,  y^ ,  can  be  obtained  when 
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the  interference  is  split  into  multiple  bands,  provided  that 
the  suppression  filter  order  is  sufficiently  large.  We  would 
expect  then  that  the  probability  of  error  performance  for 
multiple-band  interference  would  also  be  similar  to  that  given 
here  for  single-band  interference  when  the  conditions  given 
in  [2]  are  met. 

In  order  to  gain  some  insight  into  the  "bottoming  out" 
behavior  of  the  performance  obtained  using  the  Gaussian 
assumption,  we  performed  the  following  experiment.  First  we 
assumed  that  we  had  a  noiseless  estimate  of  the  transmitted 
signal  plus  interference  for  the  purpose  of  computing  the 
suppression  filter  coefficients,  but  that  the  decision  vari¬ 
able  at  the  output  of  the  sequence  correlator  is  still  cor¬ 
rupted  by  a  white  noise  term,  as  well  as  the  interference 
and  self-noise  terms.  The  resulting  error  probability  per¬ 
formance  is  the  top  curve  in  Figure  16.  We  next  assumed 
that  the  filter  coefficients  were  computed  from  a  noisy 
received  signal,  but  that  the  decision  variable  at  the  output 
of  the  sequence  correlator  was  missing  its  white  noise  term. 
The  resulting  error  probability  is  the  lower  curve  in  Figure 
16. 

We  note  first  that  the  top  curve  in  Figure  16  is  very 
little  different  from  the  performance  given  in  Figure  13  for 
the  real  world  case  where  both  the  filter  estimates  and  the 
decision  variable  are  noisy.  This  implies  that  noisy  esti¬ 
mates  of  the  interference  have  very  little  effect  on  the 
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performance  of  the  receiver  with  interference  suppression. 


At  very  low  values  of  E^/Nq  ,  the  lower  curve  falls  off 

rapidly  until  it  reaches  a  minimum  at  roughly  OdB.  The  error 

rate  increases  for  higher  values  of  E^/N^  until  it  reaches 

_  3 

a  saturation  value  of  about  3  x  10  .  This  behavior  can  be 

explained  as  follows.  At  very  low  E^/No,  the  interference 
estimate  is  very  noisy,  and  the  resulting  suppression  filter 
h&s  a  nearly  flat  frequency  characteristic,  thus  introducing 
very  little  pulse  dispersion.  As  Ej3/Nq  is  increased,  the 
interference  estimate  improves,  and  the  interference  is  sup¬ 
pressed  more  effectively.  At  about  OdB,  however,  the  self¬ 
noise  component  due  to  the  increased  dispersion  caused  by  the 
improving  interference  estimate  begins  to  wash  out  the  effect 
of  the  residual  interference.  As  the  E^/NQ  is  increased,  the 
self-noise  increases  until  the  interference  estimate  app/oaches 
its  almost  noiseless,  high  E^/NQ  value,  at  which  point  the 
error  rate  given  by  the  lower  curve  becomes  constant. 

The  conclusion  that  we  draw  from  this  experiment  is  that 
the  performance  as  given  by  the  Gaussian  approximation  is 
limited  by  the  pulse  dispersion  caused  by  the  excision  filter. 

Alternate  Bit  Error  Probability  Method  for  Symbol-by-Symbol 
Detection 

Application  of  the  Gaussian  assumption  to  the  decision 
variable,  z y  as  was  considered  in  the  previous  section,  is 
a  useful  approach  to  evaluating  the  performance  of  adaptive 
interference  cancellation  in  a  spread  spectrum  communications 
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receiver.  However,  it  is  desirable  for  several  reasons  to 
drop  this  assumption  and  investigate  other  ways  of  determining 
the  error  rate  performance  of  such  a  receiver.  The  validity 
of  the  Gaussian  assumption  is  uncertain,  and  if  we  can  find  a 
second  method  which  does  not  depend  on  this  assumption,  we  can 
get  some  idea  of  its  validity.  Moreover,  in  the  process  of 
applying  this  assumption  we  also  made  a  fairly  strong  assump¬ 
tion  about  the  statistics  of  the  PN  sequence,  i.e.,  that 

E  [c,  c.]  =0  for  k  f  j.  This  is  a  valid  assumption  for  very 
K  3 

long  sequences,  but  it  would  be  interesting  and  useful  to  be 
able  to  evaluate  the  error  probability  when  sequences  which 
do  not  satisfy  this  property  are  employed.  Finally,  the  error 
probability  expressions  given  in  the  previous  section  apply 
only  to  spread  spectrum  communications  over  nondispersive 
channels.  It  is  desirable  to  have  a  more  general  result  which 
applies  to  the  dispersive  channel  as  well. 

In  the  previous  section,  the  error  probability  was  found 
by  averaging  the  mean  and  variance  of  the  decision  variable 
over  all  possible  locations  in  the  PN  sequence,  and  assuming 
that  the  results  were  the  mean  and  variance  of  a  Gaussian 
random  variable.  In  this  section,  we  will  reverse  the  process. 
We  will  find  the  conditional  error  probability  when  a  given 
subsequence  of  the  PN  sequence  occurs  during  a  given  symbol 
transmission,  and  this  probability  wi1!  be  averaged  over  all 
possible  subsequences,  or  locations,  in  the  PN  sequence,  as 
well  as  all  possible  transmitted  sequences.  The  resulting 
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error  probability  expression  is  completely  general  in  that 
it  can  be  applied  to  any  PN  sequence  of  interest,  and  does 
not  make  any  assumptions  about  the  properties  of  the  PN 
sequence,  or  the  self-noise  caused  by  dispersion  of  the  PN 
sequence.  The  drawback  of  this  method  is  that  we  lose  the 
generality  of  the  Gaussian  assumption  and  must  now  evaluate 
the  error  rate  numerically  for  a  given  PN  code  of  interest. 

To  evaluate  the  probability  of  error  in  this  way,  we 
start  with  the  expression  for  the  decision  variable  given  . 
by  (34).  The  mean  value  of  z.  is  now  given  by: 

E(Zj)  =  2aEc  L  g(0)  u^ 
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(41) 


and  the  variance  is: 
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‘ ( zi )  =  7  E{  l  l  l  l  ck  cn  g(£)  8(m) 

^  (  k=j  L  n=jL  &=-p  m=-p 


.CO 

j  U(t)  i*  (t)  +  n(t)  n* (t) ]  p*[t  -  (k  -  £)  tc] 


p[x  -  (m  -  n)  t c ]  dt  dx 


The  effects  of  both  self-noise  ("interchip"  interference)  and 
intersymbol  interference  are  apparent  in  (41).  The  first  term 
within  the  first  brackets,  Lg(0) ,  represents  the  desired 
signal.  The  second  term  in  the  first  brackets  represents  the 
chip  pulse  dispersion  within  a  bit  interval,  the  effect  which 
we  have  been  referring  to  as  self-noise.  The  righthand  term 
in  (41)  represents  the  signal  component  due  to  chip  pulse 
dispersion  which  overlaps  the  bit  boundaries.  This  component 
is  intersymbol  interference.  We  have  made  the  assumption  in 
(41)  that  the  dispersion  caused  by  the  interference  suppression 
spans  less  than  a  single  bit  duration,  i.e.,  P  <  L. 

There  are  several  ways  of  evaluating  the  variance  of 
as  given  by  (42),  the  most  general  of  which  yields  an  ex¬ 
pression  not  much  simpler  than  (42).  Here  we  will  make  the 
simplifying  assumption  that  the  noise  and  interference  term 
of  the  sample  value  at  the  output  of  the  whitening  filter  is, 
in  fact,  white. 

In  other  words,  defining 
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P 
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and 

<J>ee  (n)  =  E[e(k)  e(k  +  n)  ]  »  (44) 

4>ee  (n)  =  °*  for  n  ^  0. 

Evaluating  4>ee(0)  yields: 
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(45) 

The  noise  and  interference  term  in  the  decision  variable 
can  now  be  written: 

(j  +  lJL-l  P 
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and  the  variance  of  is 
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Next,  define  the  normalized  noise  variance, 
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The  bit  error  probability,  conditioned  on  a  given  location  in 
the  sequence,  is 


P(error|j)  =  erfc  _1>  uj  »  uj  +  i^] 


and  the  average  bit  error  probability  is: 
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where  p„(jL  =  £)  is  the  probability  that  the  symbol  u.  starts 
at  chip  jL  in  the  PN  sequence,  and  Pu(uj.j*  uj  »  uj  +  i^  *s  t*ie 
probability  associated  with  the  subsequence  (u^^,  u^  ,  u^  +  1) 
of  the  total  transmitted  message,  u. 

Evaluation  of  (49)  depends  on  specific  knowledge  of  the 
PN  sequence  employed  and  must  be  performed  with  the  aid  of 
a  digital  computer.  Our  ability  to  do  so  is  limited  by  the 
length  of  the  sequence  in  question.  For  increasing  sequence 
lengths,  at  some  point  the  computation  time  required  to  con¬ 
sider  every  starting  point  will  become  unreasonable.  If,  on 
the  other  hand,  we  assume  that  the  sequence  is  sufficiently 
long  that  all  subsequences  of  length  L  +  2P  are  equally  likely, 
then  we  can  evaluate  (49)  by  averaging  over  all  possible 
sequences  of  length  L  +  2P.  The  processing  time  required  to 
accomplish  this,  however,  quickly  becomes  unreasonable  as 
well.  For  example,  with  a  modest  processing  gain  of  20,  and 

2  Q  g 

a  4-tap  predictor,  we  need  to  average  over  2  *  3  x  10 

sequences.  If  the  PN  sequence  in  use  is  a  maximal  length 

Y 

sequence  of  length  2  -  1,  in  order  to  satisfy  this  randomness 

property,  K  must  be  greater  than  L  +  2P. 

Equation  (49)  has  been  evaluated  for  a  maximal  length 
sequence  of  length  1023,  i.e.,  K  =  10.  Although  this  sequence 
does  not  satisfy  the  randomness  property  stated  above,  it  is 
sufficiently  long  that  the  resulting  performance  estimate  is 
quite  close  to  that  provided  by  the  Gaussian  assumption. 

Figure  17  is  a  plot  of  this  performance  estimate  for  the  same 
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interference,  suppression  filter,  and  processing  gains  as 
those  used  to  obtain  Figure  13.  Comparison  of  these  figures 
shows  that  the  performance  obtained  in  this  case  is  quite 
close  to,  but  a  bit  superior  to  that  given  by  the  Gaussian 
assumption. 

To  illustrate  the  performance  dependence  on  PN  code 
attributes,  we  have  evaluated  (49)  for  the  same  interference 
as  shown  in  Figure  17,  but  now  using  a  maximal  length  sequence 
whose  length  is  equal  to  the  processing  gain.  In  this  case 
one  complete  cycle  of  the  sequence  makes  up  the  transmitted 
waveform  for  each  bit.  Evaluation  of  (49)  was  carried  out  by 
averaging  over  each  possible  phase  shift  of  the  PN  sequence. 
This  process  was  carried  out  for  sequence  lengths  and  pro¬ 
cessing  gains  of  15,  31,  and  63;  and  the  results  are  plotted 
in  Figure  18.  It  is  evident  from  comparison  of  Figures  17  and 
18  that  use  of  the  shorter  codes  yields  superior  performance 
to  that  provided  by  the  longer  codes. 

Simulation  of  a  Bit-by-Bit  Detection  on  the  Nondispersive 
Channe l 


Another  set  of  performance  data  was  obtained  for  the 
nondispersive  channel  and  the  bit-by-bit  receiver  by  a  Monte 
Carlo  simulation.  The  conditions  which  were  simulated  were 
identical  to  those  for  which  the  calculated  performance  is 
given  in  Figure  17,  with  the  exception  that,  in  order  to 
simplify  the  simulation,  the  transmitted  information  sequence 
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in  the  simulation  consisted  of  all  "l's".  The  simulation 
processed  a  simulated  received  sequence  of  approximately 
100,000  chips,  corrupted  by  white  Gaussian  noise  and  inter¬ 
ference  consisting  of  a  sum  of  10  real  sinusoids  with  random 
phase  uniformly  distributed  over  (0,  2tt)  and  occupying  20%  of 
the  band.  The  received  sequence  was  processed  in  blocks  of 
1,200  chips;  whitening  filter  coefficients  were  computed  from 
each  block  of  1,200  samples,  using  the  least  squares  algorithm, 
and  the  1,200  samples  were  processed  using  the.  resulting 
filter  coefficients.  The  processing  of  each  block  was  over¬ 
lapped  in  such  a  way  that  there  was  no  loss  of  data  at  the 
ends  of  the  blocks,  and  phase  continuity  in  the  PN  sequence 
was  maintained  from  block  to  block  so  that  each  simulated  bit 
waveform  would  contain  as  many  subsequences  of  the  1023-chip 
PN  sequence  as  possible. 

In  order  to  obtain  a  highly  accurate  simulation,  it  is 
necessary  to  average  over  a  sufficient  number  of  samples  of 
each  possible  subsequence  of  the  PN  sequence  to  obtain  a 
good  estimate  of  the  conditional  error  probability  given  that 
the  given  subsequence  was  transmitted,  as  expressed  in  (49a). 

At  higher  SNR,  however,  it  is  usually  the  case  that  one  or  a 
few  of  the  subsequences  dominate  the  error  probability,  so 
that  it  is  sufficient  to  average  over  a  sufficient  number  of 
samples  of  each  subsequence  to  obtain  a  good  estimate  of  the 
worst  case  conditional  error  probabilities.  In  our  simulations, 
at  a  processing  gain  of  10,  we  are  processing  roughly  10,000 
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simulated  information  bits.  Since  the  PN  sequence  used  is 
1023  chips  long  there  are  1023  possible  subsequences  of 
length  10.  For  this  reason,  we  are  averaging  over  only  about 
10  samples  of  each  possible  subsequence.  Since  the  dominating 
conditional  error  probabilities  at  high  SNR  are  on  the  order 
of  10"*,  this  yields  only  a  very  rough  estimate  of  the  con¬ 
ditional  error  probabilities,  and  thus  a  rough  estimate  of  the 
total  error  probability. 

This  argument  is  borne  out  by  Figure  19,  which  is  a  plot 
of  the  results  of  this  simulation.  The  results  are  in  general 
agreement  with  those  obtained  by  the  Gaussian  approach  and  the 
conditional  probability  averaging  approach,  but  appear  more 
optimistic.  This  is  probably  a  consequence  of  the  poor 
estimates  noted  above.  A  more  accurate  simulation  could  be 
obtained  by  increasing  the  sample  size  by  perhaps  a  factor 
of  ten.  The  cost  in  computer  time  and  storage  would  be  pro¬ 
hibitive,  however,  and  it  is  not  clear  that  there  is  a  real 
need  to  do  so  since  we  have  demonstrated  a  more  efficient 
analytical  method  of  assessing  performance. 

Dispersive  Channel 

For  operation  over  a  dispersive  channel,  the  receiver 
which  performs  symbol-by- symbol  detection  will  consist  of 
the  whitening  filter  and  its  matched  filter,  as  well  as  a 
filter  matched  to  the  channel,  followed  by  the  sequence  cor¬ 
relator.  If  the  channel  is  time- varying,  the  filter  matched 
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to  the  whitening  filter,  and  the  filter  matched  to  the  channel 
can  be  realized  in  a  unit  as  the  RAKE  correlator,  as  discussed 
in  the  section  on  the  maximum  likelihood  receiver.  In  order 
to  derive  the  performance  of  this  receiver,  we  will  follow  the 
procedure  outlined  in  the  previous  section  for  the  non- 
dispersive  channel. 

As  before,  the  noise  at  the  output  of  the  whitening 
filter. 


e(k)  =  l  h (A)  .  [n(x)  +  i(r)]  p*[x  -  (k  -  JL)  x  ]  dx 

*=0  L 

(50) 


is  assumed  to  be  white.  In  this  case,  however,  the  real  part 
of  the  output  of  the  chip  pulse  matched  filter  is  not  taken, 
since  the  channel  impulse  response  is  in  general  complex, 
and  taking  the  real  part  at  this  point  would  lead  to  a  loss 
of  some  of  the  desired  signal.  The  variance  of  e(k)  is: 


<0;e(0)  -  E[|e(k)n 


P  P  f  r°° 

=  I  l  hCil)  h*(m)  [n (t)  n*  (x)  +  i(t)  i*(x)] 

1=0  m*0  JJ 


•  p*[t  -  (k  -  l)  x c ]  p [ x  -  (k  -  m)  xc]  dt  dx 


=  4Nq  Ec  g(0)  +  4EC  x c 
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The  mean  of  the  decision  variable  is 


e<v  ■  2Ec  uj  ilrxx(o)  *  2Rei  ‘X  ck  ck-»  Rxx(t)  \i 


*  2Ec  Re  Jj  Rxx“H 
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Ui  +  1  l 

3  1  k=  fi  + 1 


k-(j  +  l)L-£ 


Ck  Ck+£ 
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+  Ui-1  l  C 

3  k=j  L 


k  ckX 


where  Rxx(k)  is  as  defined  in  (25)  ,  and  L  is  the  span  of  the 
convolution  of  the  channel  with  the  whitening  filter. 

The  noise  and  interference  term  in  the  decision  variable 


‘(j+l)L-l  L 

l  l  cv  h’*(m)  e(k  +  m) 

k* j  L  m=0  K 


where  h'(m),  m  =  0,  L,  are  the  filter  coefficients  of 

the  cascade  of  the  channel  and  the  whitening  filter.  The 
variance  of  z.  is  then 


var(zp  =  \  »;e(0)[LRxx(0) 


L  (j+l)L-l 

*  2Re  J,  v  L.  Ck  Ck-”  R-(m)1  (S3) 
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We  can  form  a  signal- to-noise  ratio  from  (52)  and  (53) 
similar  to  that  in  (48)  and  use  this  expression  to  find  the 
average  bit  error  probability,  as  in  (49). 

The  resulting  expression  for  bit  error  probability 
applies  to  the  fixed  channel  case.  If  the  channel  is  time- 
varying,  then  h'(m)  and  pxx(m)  are  functions  of  time,  since 
h' (m)  is  the  convolution  of  the  channel  and  the  whitening 
filter,  and  x(t)  is  the  convolution  of  the  transmitter  pulse 
conditioning,  the  channel,  and  the  whitening  filter.  Now  the 
bit  error  probability  given  in  (49)  is  only  the  error  proba¬ 
bility  conditioned  on  a  given  channel,  i.e.. 


=  P[bit  error  |  c  (t ,  t)  ] 


(54) 


Conceptually,  at  least,  it  is  possible  to  average  the  con¬ 
ditional  probability  given  in  (54)  over  the  channel  statistics 
to  arrive  at  an  average  error  probability.  If  we  model  the 
channel  as  a  tapped  delay  line  with  time- varying  coefficients, 

M 

c  ( t ,  x )  =  l  c,(t)  6  (t  -  kx  )  -(55) 

k=l  K  c 

the  average  probability  of  error  is  given  by 


ave 


i..  j  p(c)  P[bit  errorl c (t , t) ]  dc 


(56) 


where  (56)  is  an  M- dimens ional  integral  and  p(c)  is  the 
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M-dimensional  probability  density  function  of  the  channel  tap 
weights.  Since  the  conditional  probability  has  to  be  evalu¬ 
ated  numerically,  the  integral  must  be  evaluated  numerically, 
and  in  general  the  computation  time  required  to  perform  the 
integration  would  be  completely  unreasonable. 

Using  the  technique  just  outlined  for  evaluating  bit 
error  rate  for  a  fixed  multipath  channel,  error  probabilities 
for  signaling  over  a  two-path  channel  with  tap  weights  of 
equal  magnitude  have  been  calculated  for  several  different 
interference  realizations.  In  these  calculations,  complex 
channel  coefficients  and  interference  realizations  were  used. 
Figure  20  shows  the  calculated  bit  error  rate  when  the  inter¬ 
ference  consists  of  a  sum  of  100  complex  sinusoids  occupying 
the  band  between  -0.1  and  +0.1  with  random  Hz,  phases  uniformly 
distributed  between  0  and  2tt  ,  and  the  SIR  is  -20dB.  The  PN 
sequence  used  was  the  same  maximal  length  sequence  of  length 
1023  used  previously.  The  performance  in  this  case  is  de¬ 
graded  by  roughly  lOdB  compared  to  that  obtained  with  ideal 

antipodal  binary  signaling  and  no  interference  at  an  error 

.  7 

probability  of  10  and  a  processing  gain  of  60,  and  by 
roughly  5dB  relative  to  the  nondispersive  channel  case  shown 
in  Figure  17.  Figure  21  is  the  same  channel  and  interference, 
except  that  the  SIR  is  -lOdB  here.  In  this  case  the  per¬ 
formance  is  about  2dB  better  than  when  the  SIR  is  -20dB,  when 

.  7 

the  processing  gain  is  60,  at  a  10  BER.  At  higher  values 
of  Ejj/^o  curves  *n  Figures  20  and  21  for  a  processing  gain 


Figure  20 


of  60  appear  to  converge. 

Figure  22  shows  the  calculated  performance  when  the 
channel  is  the  same  as  previously  given,  but  now  the  inter¬ 
ference  is  100  complex  sinusoids  with  uniformly  distributed 
phases  occupying  the  band  between  0.0  Hz  and  +0.1  Hz,  and  the 
SIR  is  again  -20dB.  This  figure  illustrates  that  there  is  a 
substantial  performance  improvement  when  the  interference 

bandwidth  is  reduced;  in  this  case  the  improvement  is  about 
_  7 

5dB  at  10  BER  when  the  processing  gain  is  60. 

The  interference  examples  given  here  are  near  worst  case 
for  single-band  interference  corrupting  communications  on  the 
(1,1)  channel.  This  is  because  this  channel  has  a  spectral 
peak  at  0  Hz,  and  spectral  nulls  at  ^0.5  Hz,  and  thus  the 
chosen  interference  examples  concentrate  the  interference 
power  in  the  band  where  the  signal  power  density  reaches  a 
peak.  For  this  reason,  interference  which  occupies  the  same 
bandwidth  in  other  parts  of  the  spectrum  can  be  expected  to 
cause  less  severe  performance  degradation  relative  to  ideal 
binary  communications  than  the  examples  chosen  here. 

Simulation  of  Bit-by-Bit  Detection  on  the  Dispersive  Channel 

A  Monte  Carlo  simulation  of  the  bit  by-bit  receiver  for 
the  (1,1)  channel  was  performed  in  a  manner  similar  to  that 
used  for  the  nondispersive  channel,  with  the  exception  that 
phase  continuity  in  the  PN  sequence  from  block  to  block  was 
not  maintained.  Instead,  each  block  of  1,200  samples  which 
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was  processed  in  the  simulation  started  at  the  same  location 
in  the  PN  sequence,  and  as  a  result  only  a  subset  of  all 
possible  subsequences  in  the  PN  sequence  were  transmitted  in 
the  simulation.  We  would  expect  that  for  this  reason  the 
complete  set  of  worst  case  subsequences  was  not  used,  and 
the  resulting  performance  estimate  would  be  optimistic. 

Figure  23  shows  the  results  of  this  simulation.  The  inter¬ 
ference  used  in  this  simulation  was  a  sum  of  10  complex 
sinusoids  with  random  phases  occupying  101  of  the  band  (0.0  Hz 
to  0.1).  Comparing  Figure  23  with  Figure  22  confirms  that 
the  simulated  performance  exhibits  generally  the  same  sort  of 
behavior  as  the  calculated  performance,  although  somewhat  more 
optimistic. 

A  Monte  Carlo  simulation  of  the  bit-by-bit  receiver  for 
fading  multipath  was  also  performed.  This  receiver  was 
realized  using  a  RAKE  correlator,  and  simulations  have  been 
performed  using  both  an  adaptive  RAKE  and  a  fixed  RAKE.  The 
adaptive  RAKE  simulation  was  implemented  using  a  "snapshot" 
approach.  In  this  technique  a  set  of  channel  coefficients 
are  chosen  pseudo- randomly  and  a  block  of  2,000  chips  are 
transmitted  through  the  resulting  channel.  At  the  receiver, 
the  interference  suppression  filter  coefficients  are  computed 
from  the  block  of  data.  The  first  1,000  chips  are  then  used 
to  train  the  RAKE  and  the  next  1,000  chips  are  used  for  col¬ 
lecting  error  rate  data.  The  simulation  of  the  fixed  RAKE 
was  performed  in  a  similar  manner  except  that  the  full  2,000 
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chips  in  each  block  were  used  for  collecting  error  rate  data, 
and  the  exact  matched  filter  coefficients  were  used  for  the 
RAKE  coefficients  instead  of  estimates  obtained  from  the  RAKE 
adaptation  process.  This  method  has  the  advantage  of  providing 
a  larger  sample  size  for  the  purpose  of  estimating  the  error 
rate.  Both  simulation  methods  were  used,  and  it  was  found 
that  the  adaptive  RAKE  exhibited  performance  which  was  only 
slightly  degraded  from  that  of  the  fixed  RAKE. 

Figure  24  shows  the  results  of  the  fixed  RAKE  simulation, 
using  100,000  chips  of  data  processed  in  50  blocks  of  2,000 
chips  each.  Thus  there  were  50  channel  realizations.  The 
channels  had  two  active  taps,  the  first  and  the  fourth,  which 
were  chosen  pseudo- randomly  from  a  complex  Gaussian  distri¬ 
bution  with  zero  mean  and  unit  variance.  The  interference 
was  a  sum  of  20  complex  sinusoids  with  random  phases  occupying 
20%  of  the  band  (-0.1  to  +0.1  Hz).  Also  plotted  in  Figure  24 
is  performance  of  ideal  dual  diversity  communications  with  no 
interference.  It  is  seen  from  this  figure  that  the  perform¬ 
ance  in  the  presence  of  interference  tracks  that  of  ideal 
dual  diversity,  except  that  it  is  degraded  by  a  fjw  dB. 

Figure  25  shows  the  simulated  performance  for  the  same  set 
of  channels,  but  with  interference  which  covers  101  of  the 
band  (0.0  to  0.1  Hz).  The  performance  in  this  case  is 
improved  over  that  of  the  201  interference,  as  would  be  ex¬ 
pected. 

All  the  limitations  of  our  simulation  technique  which 
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were  stated  for  the  nondispersi ve  channel  and  the  fixed 
multipath  channel  apply  even  more  strongly  for  this  simula¬ 
tion  since  now  we  need  not  only  average  over  the  noise  and 
interference  density  functions  and  the  subsequences  of  the  PN 
sequence,  but  we  must  average  over  the  complex,  two-dimensional 
probability  distribution  of  the  channel  as  well.  Since  we  have 
not  increased  the  sample  size  over  that  of  the  fixed  channel 
simulation,  we  would  expect  that  the  performance  estimate 
given  by  the  simulation  would  be  less  accurate  than  that  given 
by  the  fixed  channel  simulations.  For  this. reason  it  is  not 
surprising  that  a  few  of  the  data  points  in  these  simulations 
fall  below  the  double  diversity  performance.  This  is  pre¬ 
sumably  a  consequence  of  the  limitations  of  the  simulation. 

A  considerable  increase  in  sample  size  would  be  required  to 
provide  an  accurate  estimate  of  performance,  and  unfortunately, 
the  resulting  storage  and  computation  time  requirements  would 
be  prohibitive.  However,  the  simulations  which  we  have  per¬ 
formed  indicate  that  adaptive  interference  suppression  tech¬ 
niques  can  be  used  with  good  results  in  a  fading  multipath 
environment . 

IV.  Error  Probability  Bounds  for  Maximum  Likelihood  Detection 

In  order  to  establish  error  rate  bounds  for  maximum 
likelihood  sequence  detection  of  direct  sequence  spread 
spectrum  signals  with  receiver  processing  of  the  type  con¬ 
sidered  in  this  report,  we  first  consider  the  probability  of 
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a  pair-wise  sequence  error.  This  is  the  probability  that  the 
transmitted  information  sequence,  u,  is  decoded  at  the  receiver 
as  the  sequence  u' ,  where  u'  is  identical  to  u  in  all  but  a 
given  number  of  bit  positions.  We  then  make  use  of  this 
pair-wise  sequence  error  probability  to  establish  bounds  on 
the  bit  error  probability.  The  arguments  given  here  follow 
fairly  closely  those  given  in  [5]  for  probability  of  error 
on  an  intersymbol  interference  channel. 

In  order  to  find  the  probability  of  a  pair-wise  sequence 
error,  we  start  by  defining  the  decision  variable  z  =  A  -  A', 
where  A  is  the  likelihood  function  given  in  (19)  for  the 
correct  transmitted  sequence  decision  u  and  A'  is  the  likeli¬ 
hood  function  for  some  incorrect  sequence  decision,  u'.  Thus 
the  probability  of  a  pair-wise  error  is  given  by 

Pr(£)  =  Pr(z  <  0) 

where  e  is  the  error  sequence  e  *  (u  -  u').  Again  we 
assume  that  the  decision  variable  z  is  a  Gaussian  random 
variable,  so  that  we  only  need  to  find  the  mean  and  variance 
of  z  to  have  the  desired  probability  of  error.  The  mean  of 
z  is 
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Substituting  the  error  sequence  c  into  (57): 
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We  will  next  evaluate  the  variance  of  z  on  the  assumption, 
used  in  the  previous  section,  that  the  total  noise  at  the 
whitening  filter  output,  due  to  the  white  noise  and  the  narrow- 
band  interference,  is  white.  Thus  the  noise  term  in  z  is 
given  by: 


N 


7 


NL-1  P 

2  I  I  C u  I  k  /  l! 

k=-NL  £  =  0  LK/LJ 


U'[k/Lj)  Ck  h(i)  e(k  -  ^ 


(59) 


and  the  variance  of  z  is  given  by 


var(z)  =  E(NZ) 


NL-1  P  NL-1  P 


4  J«L  Jo  n=-NL  Jo  tULk/LJ  '  U'^ 


Cu^n/LJ  “ u*  Ln/Lj)  ckcnh(A)  h(m)  Ete(k"^  e(n-m)] 
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NL-1  P  P 

'  16  k  =  ?NI.  Jo  Jo  ELk/LJ  C  (_(k+m*»)/Lj 

'  ck  ck*m-l  h(1)  h(n,)  »ee(0) 

N-  1 

-  I  [LCe,)2  g(O) 

j  =  -  N  J 

P  ( j  + 1) L- 1 

'  \h  JjL  Ei  «k  «*-.  «W1 

(60) 

Finally,  defining  the  signal- to-noise  ratio  as 
y  =  [E(z)]2/2  var(z) 

we  have 

a2  E  N-l  - 

Y  =  - l  [L( I-V  g(0) 

ae  j-N 

P  (j+lJL-l 

+  2  m=l  k  =  j  L  ^  £L(k-m)/Lj  Ck  ck-m  «(m)]  C61) 

2 

where  a ‘  is  given  by  (47).  Thus  the  probability  t’  it  the 
maximum  likelihood  receiver  will  select  the  sequence  u’  instead 
of  u  is  given  by: 
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P(e)  =  \  erfc[/Y] 

In  order  to  compute  the  bit  error  probability  we  start 
by  defining  the  error  sequence  £^  as  any  error  sequence  which 

■f  Vi 

has  a  non-zero  element  for  the  first  time  in  the  *  position, 
i.e.,  =  0  for  j  <  i .  Thus  we  can  union  bound  the  bit  error 

probability  given  that  the  first  error  occurs  at  the  Jl  node 
by : 


V*)  £  l  2  P(£A)  (62) 

The  summation  in  (62)  is  over  all  possible  error  sequences 
which  begin  at  node  i ,  and  w(e^)  is  the  number  of  places  in 
which  is  non- zero.  In  order  to  obtain  a  union  bound  on 
the  total  probability  of  bit  error,  we  also  must  average  over 
all  nodes  in  the  sequence.  This  is  because  for  a  given  non¬ 
zero  error  sequence  starting  at  two  different  nodes,  the  value 
of  P(£)  is  different,  depending  on  the  particular  subsequence 
of  the  PN  sequence  which  occurs  at  that  node.  For  example, 
the  probability  that  a  single  bit  error  will  occur  at  node  p 
depends  on  the  PN  sequence  on  the  interval  pL  £  k  <  (p  +  1)L, 
while  the  probability  that  a  single  bit  error  will  occur  at 
node  q  depends  on  the  PN  sequence  on  the  interval  qL  £  k  < 

(q  +  1)L.  In  general,  the  PN  sequence  in  these  two  intervals 
is  not  the  same  unless  the  processing  gain,  L,  is  equal  to 
or  some  integer  multiple  of  the  PN  sequence  length. 

J 
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The  bit  error  probability  is  then  bounded  by: 


Pb  1  l  l  w(££)  2  W(i£)  p(££)  =  l  P(£)  PbU)  (63) 

e£  i 

The  outer  summation  is  over  all  possible  nodes  at  which  an 
error  event  can  start,  or  equivalently,  over  all  points  in 
the  PN  sequence  at  whu  the  bit  u^  can  start.  The  proba¬ 
bility  p (£)  is  the  proL  bility  that  u^  starts  at  a  particular 
point  in  the  PN  sequence.  We  assume  that  this  probability 
is  uniformly  distributed  over  the  PN  sequence,  or  p ( it )  = 

1/ (sequence  length). 

We  will  start  by  evaluating  the  conditional  probability 
Pfo(l).  To  do  so,  we  will  need  the  inequality:  erfcVx  £ 
erfcVy{exp(y  -  x}}  [6].  We  will  also  define  y^fe^)  as  the 
signal- to-noise  ratio,  y,  given  by  (61),  evaluated  for 
e ■  -  0 ,  j  <  i ,  and  e  =  +  1.  We  make  the  further  definition 

J  ** 

that 


a2  E 


Jl,min 


_  P  (H+l)m- 1 

-  [Lg  (0)  *  2  l  l  c,  c.  g(m)  ] 


m=l  k=£L+m 


By  Virtue  of  these  definitions  and  the  above  inequality,  we 
can  write  the  following  inequality: 


must  be  evaluated 


In  order  to  evaluate  (65),  the  term  e 
using  an  error  state  diagram  and  solving  for  its  transfer 
function.  This  technique  was  developed  by  Forney,  Viterbi  and 
Omura,  and  others  [5,6],  for  the  purpose  of  finding  error 
probabilities  for  intersymbol  interference  channels.  Ex¬ 
tending  this  technique  to  (65)  is  a  bit  problematic,  however, 
since  y ^  depends  not  only  on  the  given  error  sequence,  but 
also  on  the  timing  of  the  PN  sequence  relative  to  the  error 
sequence.  Another  way  of  stating  the  problem  is  that  although 
we  are  employing  a  binary  antipodal  signaling  scheme,  the 
basic  channel  waveform  varies  from  bit  to  bit.  For  this 
reason,  the  resulting  error  state  diagram  will  have  time- 
varying  path  gains.  Although  we  know  of  no  exact  solution 
to  such  a  time-varying  flowgraph,  we  wil]  proceed  to  construct 
the  flowgraph,  and  write  an  approximate  solution. 

First  we  will  substitute  (65)  into  (62)  : 

PK(£)  <  \  erfc/y„  _•  exp(Y.  .  )  7  w(e„) 

b v  J  —  2  '£,  min  K'-'£,mm'  L  K—Z 

-Z 

■  2  W(£jl)  expt-y^ep] 


=  7  erfc/Y*,min  exP <YliBin)  l 

— £ 

N-l  ,  a2  E 

•  TrfVrr  exP[ - r5  E(£n-»)  (««) 

J-£-  2  J  ae 
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where 


E(e£>j)  =  L (e j ) 2  g(0) 

P  Cj+1)L-1 

+  2  nJl  k=ZjL  ^  £L(k-m)/Lj  ck  Ck-m 

(67) 

We  have  assumed  that  an  error  occurs  for  the  first  time  at 
node  Z,  i.e.,  e.  =  0,  j  <  Z,  e0  =  + 1 .  Thus,  for  j  <  Z 
E(e£,j)  =  °>  and 

P  (A+1JL-1 

E(e,,  *)  =  Lg(0)  ♦  2  l  I  ck  ck  g(m) 

m-1  k=£L+m  K  K  m 

(68) 


and  we  define  the  quantity  aQ (A)  as 


a0(A)  =  exp  [- 


a 


E(e£,0] 


(69) 


There  are  three  possible  values  for  E(£^,  £  +  i) ,  for  i  £  1  and 
e j^O ,  Z  <  j  <  Z  +  i.  If  eJl+i  =  e£+i-1,  then: 

’  P  (A+i+1)  L- 1 

E(e  A+i)  =  Lg(0)  +  2  \  l  ck  c.  g(m) 

m-1  k- (A+i)L  K  K  m 


(70) 
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and  we  define  the  quantity  a^(H  +  i)  by 

2  p 

a  Ec 

a,  (A  +  1)  =  exp  [-  - 5—  E(e  ,  £+i)] 

'e 


If  eJt+i  "e£+i-l* 


then 


Etc.,  i*i)  -  Lg(0)  *  2  l  [  W  + 

"-1  k=(i+l)L+m  k 


(je,+i)L+m- 1 

I  cv  ck-m]  gO) 

k= (H+i) L  K  K  m 


and  we  define  the  quantity  a2(£  +  i)  by: 


a2  E 

a2(Ji  +  i)  =  exp  [-  - -  Efe^,  *+i)] 

a 


If  e£+i  =  °»  then  A  +  i)  =  °* 

The  resulting  flowgraph  is  shown  in  Figure  26.  The 

1  -w(e.) 

factors  of  in  each  branch  account  for  the  factors  2  J 

in  (66),  and  the  factor  1  is  also  inserted  into  each  branch 

to  indicate  that  a  bit  error  is  made  on  every  transition  in 

the  error  state  diagram.  The  variable  i  is  incremented  each 

time  a  non- zero  state  is  passed  through  on  a  y-path  from  the 

start  node  to  the  finish  node.  is  obtained  by 
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differentiating  the  resulting  generating  function  with 
respect  to  I,  as  shown  by  Forney,  and  Viterbi  and  Omura  [5,6]. 

As  was  previously  mentioned,  it  is  not  possible  to  write 
a  closed  form  expression  for  the  transfer  function  of  this 
error  state  diagram,  since  the  path  gains  are  time-varying. 
However,  because  the  transfer  function  is  the  sum  of  all 
possible  paths  through  the  state  diagram,  we  can  write  the 
first  several  terms  of  the  transfer  function  as  follows: 

T(a0,  a1,  a2,  l ,  I)  =  aQ(£)I  +  a^fc  ♦  1)  aQ(lt,)  1 2/ 2 

♦  a0(£)  a2U  +  1)  1 2 / 2 

+  a0(*)  a^l  *  1)  a^Jt  +  2)  I3/4 

♦  a0(£)  ax(£  +  1)  a2(Jt  +  2)  1 3/ 4 

+  aQ  (£.)  a2(l  +  1)  3l(£  +  2)  I3/4  +  ...  (74) 


We  have  written  in  (74)  terms  representing  all  possible  ways 
of  making  one,  two,  and  three  bit  errors.  In  order  to  find 
an  approximate  upper  bound  for  the  bit  error  rate,  we  take 
the  derivative  of  (74)  with  respect  to  I  and  set  1*1. 


«  aQ(A)  +  aQ(£)  a^ (£  +  1) 


+  a 0(£)  a2(*  +  1)  ♦  |  aQ(£)  ai^£  +  D  ai^£  +  2)  + 


(75) 
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Replacing  (75)  in  (66),  we  have: 

Pk(2)  <  4  erfc/y*  ~  exp(v.  .  )  [a  (2.) 
bl  y  -  2  r2,min  hUH,niinn  1 

+  aQ(2)  ax(2  +  1)  +  aQ (2)  a2(2  +  1) 

+  |  aQ(2)  ai(2  +  1)  ax(2  +  2)  +  ...]  (76) 

a2  Ec 

Noting  that  Y£>m.n  =  -2-  E(£j1,  2),  and  thus  exp(Yjl>min)  = 
[ag (2) ]  1 ,  we  have  “ 

V*>  <  |  erfc/i^—  [1  ♦  aj(*  ♦  1)  .  a .,(«  <■  1) 

*  |  ajCH  *  1)  EjCt  *  2)  ♦  ...]  (77) 

Dropping  all  but  the  first  term  in  (77)  gives  the  approximate 
upper  bound: 


Pb  (2)  <  -j  erfc/y 


2  ,min 


(78) 


This  represents  the  probability  of  error  due  to  the  minimum 
distance  path  through  the  error  state  diagram,  assuming  that 
the  error  event  starts  at  node  2. 

Finally,  to  find  the  average  bit  error  probability,  we 


substitute  (78)  into  (63),  yielding: 
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Pb  <  1  I  pU)  erfc/v7—  <79> 

This  expression  has  been  evaluated  for  the  nondispersive 
channel  and  several  PN  sequences.  Figure  27  is  the  resulting 
calculated  performance  of  the  maximum  likelihood  receiver 
when  the  interference  is  a  sum  of  100  real  sinusoids  occupying 
the  band  between  0.0  and  0.1  Hz,  with  independent  and  identi¬ 
cally,  uniformly  distributed  phases  between  0  and  2 it .  The 
SIR  is  -20dB,  so  the  performance  shown  is  for  the  same  channel 
and  interference  as  is  illustrated  in  Figure  17  for  the 
receiver  which  performs  bit  decisions  on  a  bit-by-bit  basis. 
Comparison  of  Figure  27  with  Figure  17  reveals  that  the 
receiver  which  performs  maximum  likelihood  sequence  estimation 
has  only  slightly  improved  performance  over  the  bit-by-bit 
receiver,  particularly  for  the  high  processing  gain.  This  is 
true  because  at  the  higher  processing  gain,  the  inequality 
L  >>  P  holds,  i.e.,  the  symbol  duration  is  much  greater  than 
the  dispersion  caused  by  the  interference  suppression  filter, 
and  for  this  reason  there  is  negligible  intersymbol  inter¬ 
ference. 

The  remaining  performance  degradation  relative  to  ideal 

.  7 

binary  signaling  without  interference,  about  5dB  at  10  BER, 
is  due  to  residual  interference,  loss  of  signal  power  due  to 
the  notch  created  by  the  interference  suppression  filter,  and 
the  chip  pulse  dispersion  or  interchip  interference.  The 
loss  of  signal  power  due  to  the  suppression  filter  is  20%,  or 
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ldB,  since  the  interference  occupies  20%  of  the  band. 

Assuming  that  the  residual  interference  at  the  output  of  the 
matched  filter  is  negligible  relative  to  the  white  noise, 
the  remaining  4dB  of  degradation  at  10  BER  must  be  due  to 
interchip  interference.  Furthermore,  since  the  performance 
curve  for  a  processing  gain  of  60  continues  to  diverge  from 

_  7 

the  ideal  binary  curve  below  10  BER,  it  appears  that  the 
interchip  interference  increases  with  increasing  E^/Nq. 

Figure  28  shows  the  performance  with  the  same  channel 
and  interference,  but  with  an  SIR  of  -lOdB.  Here,  at  the 
higher  processing  gain,  the  performance  is  degraded  by  only 
1.5- 2dB  at  10  7  BER.  Thus  at  lower  levels  of  interference, 
it  appears  that  the  suppression  filter  contributes  signifi¬ 
cantly  less  interchip  interference  than  -20dB,  in  this  case 
causing  less  than  ldB  degradation  above  that  caused  by  loss 
of  signal  power. 

The  performance  bounds  given  in  Figures  27  and  28  were 
obtained  with  a  maximal  length  sequence  of  length  1023.  We 
next  consider  the  performance  with  shorter  sequences.  Figure 
29  is  the  performance  obtained  when  the  channel  and  inter¬ 
ference  are  the  same  as  given  in  Figure  27,  but  the  sequence 

length  is  now  63.  We  see  that  at  the  highest  processing  gain, 

- 1 2 

the  performance  at  error  rates  above  10  is  not  much  dif¬ 
ferent  than  that  obtained  with  the  longer  sequence;  however, 
at  the  lower  processing  gain,  the  BER  bottoms  out  at  a  lower 
error  rate  when  the  shorter  sequence  is  used. 
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Finally,  Figure  30  shows  the  performance  obtained  when 
the  sequence  length  equals  the  processing  gain,  for  sequence 


lengths  of  15,  31,  and  63.  Here  again,  at  the  high  processing 
gain,  the  performance  is  not  much  different  from  that  obtained 
with  the  long  spreading  sequence,  but  at  lower  processing 
gains,  the  minimum  achievable  error  rates  are  significantly 
better  than  those  obtained  with  the  long  PN  sequence. 

Dispersive  Channel 

The  upper  bound  on  error  rate  performance  of  the  maximum 
likelihood  receiver  for  transmission  over  a  nondispersive 
channel,  which  is  summarized  in  (65)  and  (79)  can  be  modified 
to  give  an  upper  bound  on  performance  of  the  maximum  likeli¬ 
hood  receiver  for  transmission  over  a  dispersive  channel. 

This  is  done  in  a  manner  similar  to  that  used  to  generalize 
the  alternate  method  for  evaluating  the  bit  error  rate  per¬ 
formance  of  the  bit-by-bit  receiver  given  in  Part  III.  This 
has  been  done,  and  the  results  have  been  used  to  evaluate  the 
upper  bound  for  channels  with  weights  (1,1),  (.5,. 71,. 5)  and 
(. 38, .6, . 6, . 38) .  These  are  channels  given  by  Magee  and 
Proakis  [7]  as  the  channels  with  length  2,  3,  and  4  which 
have  minimum  distance  of  all  channels  of  their  respective 
lengths  for  binary  antipodal  signaling  in  white  Gaussian 
noise.  Although  these  channels  will  not  necessarily  yield 
minimum  distance  in  non-white  noise,  they  are  useful  as 
examples  of  the  performance  of  the  maximum  likelihood  receiver 
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for  channels  with  severe  multipath  distortion. 

Figures  31  and  32  show  the  performance  of  the  maximum 
likelihood  receiver  for  the  channel  with  weights  (1,1)  and 
interference  which  consists  of  100  complex  sinusoids  with 
random  phases  uniformly  distributed  between  0  and  2ir ,  and 
occupying  the  band  between  -0.1  and  +0.1  Hz.  The  PN  sequence 
used  in  these  examples  is  the  maximal  length  sequence  of 
length  1023  used  in  previous  examples.  Figure  31  illustrates 
the  performance  with  -20dB  SIR,  32  with  -lOdB  SIR.  Comparing 
these  results  with  the  corresponding  results  for  bit-by-bit 
detection  (Figures  20  and  21)  shows  that,  as  with  the  non- 
dispersive  channel,  at  a  processing  gain  of  60,  there  is  very 
little  difference  in  the  performance  of  the  maximum  likeli¬ 
hood  and  bit-by-bit  receivers.  At  the  lower  processing  gains 
the  maximum  likelihood  receiver  again  bottoms  out  at  a  lower 
error  probability  than  the  bit-by-bit  receiver. 

The  performance  of  the  maximum  likelihood  receiver  for 
the  (1,1)  channel  when  shorter  sequences  are  used  was  also 
considered,  and  the  results  were  similar  to  those  *  ined 
for  shorter  sequences  with  a  nondispersive  channel.  a 

processing  gain  of  60,  the  performance  differs  little  from 
that  of  the  longer  sequence.  At  smaller  processing  gains, 
the  bottoming  out  behavior  is  considerably  less  severe  than 
that  exhibited  when  the  longer  (1023  chips)  sequence  is  used. 

Figure  33  shows  the  performance  for  the  same  channel,  inter- 

» 

ference,  and  processing  gains  as  Figure  31,  but  with  a  63-chip 
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PN  sequence.  Figure  34  is  the  performance  with  same  channel 
and  interference,  but  now  the  processing  gain  and  sequence 
length  are  the  same,  with  values  of  15,  31,  and  63. 

Finally,  Figures  35  and  36  show  the  performance  obtained 
with  the  3-  and  4- tap  channels  given  above.  In  these  compu¬ 
tations  the  interference,  sequence  length  and  processing 
gains  are  the  same  as  those  which  yielded  the  performance 
curves  in  Figure  31.  These  curves  exhibit  behavior  similar 
to  that  of  the  nondispersive  and  (1,1)  channel  except  that  as 
the  span  of  the  dispersion  increases,  the  degradation  relative 
to  ideal  binary  antipodal  signaling  increases,  reaching 

-  7 

approximately  16dB  at  10  BER  for  the  4- tap  channel  with  a 
processing  gain  of  60.  Although  this  is  a  severe  degradation, 
we  have  demonstrated  that  the  performance  of  the  maximum 
likelihood  receiver  using  autoregressive  spectral  estimation 
techniques  for  suppressing  narrowband  interference  in  an 
environment  of  severe  multipath  distortion  exhibits  error 
probability  performance  which  decreases  exponentially  as  a 
function  of  Eb/NQ.  This  is  a  significant  improvement  over 
performance  which  is  achievable  even  on  a  nondispersive 
channel  without  interference  suppression. 
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>endix  A 


In  a  previous  report  [2],  the  signal- to-noise  ratio, 
2 

defined  as  E  (z)/var(z),  of  the  decision  variable  at  the 
output  of  the  sequence  correlator  shown  in  Figure  11  is 
shown  to  be: 


SNR0  =  p  p  p  P 

l  h2(t)  +  l  l  h(£)  h(m)  p.  (£-m)+o2  £  h2(£) 

£=1  1=0  m=0  £  =  0 

(A.  1) 

This  expression  is  limited  for  the  following  reasons:  (1)  it 
only  applies  when  the  suppression  filter  does  not  delay  the 
desired  signal,  i.e.,  the  first  tap  is  the  largest  tap;  (2) 
it  only  applies  to  baseband  signaling.  In  this  appendix  we 
derive  a  more  general  formula  for  the  signal- to-noise  ratio 
for  bandpass  signaling,  with  a  zero  delay  suppression  filter 
and  with  a  linear  phase  suppression  filter. 

When  the  whitening  filter  alone  is  used,  the  suppression 
filter  has  zero  delay,  and  the  decision  variable  is  given  by 
(35): 

(j  +  lJL-i  P 

zj  •  2“Ec  Lh'°’  uj  *  2oEc  uL(1c-1)/lJ 


•  ck  ck-t  hW 
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(j+lH-i  P 

*  I  I  clf  h(£)  Re 

k  =  j  L  £  =  0  K 


[n(x)  +  i  (t)  ] 


p*[x  -  (k  -  £)  tc]  dx 


(34) 


The  mean  of  z ^  is 


E(z .)  =  2aEc  Lh (0) 


(A.  2) 


and  the  variance  is 


var(Zj)  -  S.  *  Ij  *  Nj 


(A. 3) 


The  first  term  in  (A. 3)  is: 


2  2  U+1JL-1  P  P 

S,  =  4aZ  Ej  l  I  ll 

J  k= j  L  m= j  L  £  =  1  n=l 


’  u\(k-£)/Lj  u|_(m-n)/Lj  ck  cm  ck-£  Cm-n  h(£)  h(m) 

(A. 4) 

Since  {c^}  is  assumed  to  be  a  white  sequence,  c^  and  c^  are 

uncorrelated  unless  k  =  £.  Thus  c^  c^  cm  cn  =  0  unless  k  *  £ 

and  m  *  n,  k  ■  n  and  £  =  m,  or  k  =  m  and  £  =  n,  in  which 

case  c7~cT~c~~c~  =  1.  So  c,  c  c,  0  c  =  0  unless  k  =  m 
k  £  m  n  km  k-£  m-n 

and  £  *  n,  k  *  m  -  n  and  m  =  k  -  £,  or  k  =  k  -  £  and 
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m  =  m  -  n.  The  first  condition  is  easily  satisfied.  The 
second  condition  is  two  equations  in  three  unknowns  which 
yields  the  following  equality: 

1  -  -n 


which  cannot  be  satisfied,  since  the  intervals  of  summation 
of  both  1  and  n  are  positive.  Using  the  first  condition  in 
(A. 4) ,  we  get 


.  =  4a2  L  l  h‘(l) 
3  c  1=1 


(A.  5) 


The  second  term  in  (A. 3)  is: 


(j  +  lJL-1 


«*VL'1  ?  ?  _ 

l  l  l  ck  cm 

m= j  L  1=0  n=0  K  m 


J{  i (t)  r*TTJ  p*[t  -  (k  -  1)  Tc] 


*  p[t  -  (m  -  n)  t  ]  dt  dr 


1  ?  f  (j-UL-1 

7  11 

L  1=0  n=0  k 


r  If 


-  T) 


•  P*[t  -  (k  -  1)  xc)  p[T  -  (k  -  n)  Tc]  dt  dt 


(A. 6) 
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Equation  (A. 6)  is  evaluated  in  a  manner  similar  to  that  used 
to  evaluate  (28),  yielding 


P  P 


I.  -  2L  x  E  Y  I  h(£)  h(n)  <f>-  •  [  (n  -  t)  x  ] 
3  c  c  £=0  n-0  11  c 


(A. 7) 


The  third  term,  ,  can  be  found  in  a  straightforward  manner 
to  be: 


N.  =  2L  N  E  l  h2(Z) 

J  o  C  J  n 


(A. 8) 


Summing  (A. 5),  (A. 7),  and  (A. 8)  gives  the  expression  for 
var(Zj)  given  by  (37): 


p  ? 

var(z-)  =  2L  N  E_  [  h 
3  °  c  £=0 


P  P 

+  2L  x  E  l  l  h  (£)  h  (m)  <j>..[(m  -  l)  x  ] 
c  £®0  m=0  11  c 


♦  4a2  LE2  l  h2  (£) 
C  £=1 


The  resulting  SNR,  now  defined  as  E  (z.)/2  var(z.)»  is: 

J  3 


L  Ec  h2 (0)  a2 

- p — - p — p - - - p — ; 

N  l  h2 (£)+x  l  l  h(£)h(m)<j>.  .  [(m-fc)x  ]  +  2aZE  [  h' 

0  £*0  £*0  m»0  li  c  c  £-1 
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This  expression  is  similar  to  (A.l)  with  the  exception  that 
(38)  is  a  bit  more  general  than  (A.l),  and  that  when  p^fO)  * 
*ii(0),  the  interference  term  in  (38)  contributes  relatively 
half  as  much  to  the  denominator  in  (38)  as  it  does  in  (A.l). 
This  is  because  a  coherent  bandpass  receiver  has  a  3dB  ad¬ 
vantage  over  a  baseband  receiver  against  interference. 

Next  we  consider  the  linear  phase  filter  case.  For  con¬ 
venience,  we  allow  the  linear  phase  filter  to  have  an  acausal 
impulse  response.  The  linear  phase  interference  suppression 
filter  has  coefficients  g(t),  i  =  -P,  ...»  P,  where  g(£)  = 
g(-£),  and  the  g(£)'s  are  real  valued.  The  decision  variable 
in  this  case  is  given  by: 


(j  +  llL-1  P 

z.  .  2a  Ec  Lg CO)  Uj  *  2a  Bc  J 

l?0 

•  Ck  Ck-£  *<*) 


+ 


(j  +  l)L-l 
k = j  L 


8(A) 


r 
-  00 


[n(T)  +  i(x)  ] 


•  p*[x  -  (k  -  i )  t c  ]  dx 


(34) 


The  mean  of  this  decision  variable  is 


E(z.)  -  2a  Ec  Lg(0) 


(A.  9) 


103 


and  the  variance  is  again  given  by 

var(Zj)  =  Sj  +  1^  +  Nj  (A. 10) 

The  second  two  terms  in  CA.10)  are  found  in  a  manner  similar 
to  that  used  to  find  the  second  two  terms  in  (A. 3).  These 
terms  are: 


P  P 

I.  *  2L  r  E  l  l  gtfc)  g(n)  <J>.  [ (n  -  l)  t  ] 

J  l*-p  n=-p  11  c 

(A. 11) 


and 


N  =  2L  N  E  l  g2(t) 
J  °  c  t=-p 


(A. 12) 


The  self-noise  variance,  ,  is  found  as  follows: 


Sj  •  4a2  e| 


(  j  + 1 )  L  -  1  P  (j  +  l)L-  1  P 

I  l  1  l 

k= j  L  Jl=-p  m=  j  L  n=-p 

ij^O  n^O 


’  u|_(k-t)/LJ  u^(m-n)/Lj  ck  ck-t  cm  cm-n  gU)  g(n) 

(A. 13) 


The  three  conditions  outlined  above  under  which  c,  c.  0  c  c 
■  0  apply  equally  well  for  (A. 13).  However,  in  this  case,  the 
first  condition  as  well  as  the  second  condition  apply,  i.e., 


we  can  set  k  »  m  and  i  »  n,  but  we  can  also  set  m  ■  k  -  l  and 
n  *  -l.  Thus  we  can  write: 


,  2  (J*1JL-1  P  2 

i  K  I  l  sW 

k=jL  £*-p 
t/0 


2  2  ( j  +  1J  L- 1  P 

+  4a2  E2  I  l  g(£)  gC-A) 

k= j  L  £=-p 


jL<k-  £<(j  +  l)L-l 


(  P  P 

.2  c2  )  r  _2/,i  .  ->  v 


4a2  E‘  2L  l  g 2(t)  +  2  £  (L  -  £)  g2(i) 


8a2  E2  I  (2L  -  l)  g2U) 
c  1=1 


(A. 14) 


Summing  these  three  terms  yields  the  expression  for  var(z^) 

2 

given  by  (39b),  and  the  ratio  E  (Zj)/2  var(Zj)  is  given  by 
(40). 

An  SNR  improvement  factor,  defined  as 


SNR  with  filter 


was  given  in  [2]  as 


10S 


T - F 


Pi(°)  +  on 


l  hz(A)+  I  I  h(t)  h (m)  p.  (£-m)  +  l  h£(i) 
L=1  1=0  m=0  1  n  fc=0 


(A. 15) 


A  similar  expression  for  improvement  factor  can  be  obtained 
by  taking  the  ratio  of  (38)  over  (33): 


[No  +  Tc  ♦iiC0)]  h2(°) 

n  *  p  p  p  p 

N  I  h2(t)  +  T  l  l  h(£)h(m)*..  [(m-£)tr]  +  2aZE2  f  h2(t) 

0  i=0  C  1=0  m=0  11  c  c*=l 

(A. 16) 

If  h(k) ,  k  *  0,  . . . ,  P  is  normalized  so  that  h(0)  *  1,  and  we 
let  on  =  NQ/2»  a  =  1,  tc  =  1,  and  Ec  =  1,  then  (A. 16)  reduces 
to  (A. 15),  with  the  exception  that  the  interference  terms  in 
the  numerator  and  denominator  are  a  factor  of  2  smaller  than 
in  (A. 15).  This  is  due  to  the  fact  that  (A. 15)  applies  to 
the  baseband  case,  and  (A. 16)  to  bandpass  signaling.  The 
coherent  bandpass  receiver  thus  has  a  3dB  advantage  over  the 
baseband  receiver  against  interference,  as  has  been  pointed 
out  earlier. 

The  improvement  factor  given  by  (A. 16)  applies  to  the 
receiver  with  a  zero  delay  interference  suppression  filter. 

The  improvement  factor  for  the  receiver  which  contains  a 
linear  phase  suppression  filter  can  be  obtained  by  taking  the 
ratio  of  (40)  over  (33) : 


g(A)g(m)^ii[Cm-i)Tc 

p 
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Appendix  B 


In  view  of  the  results  of  Appendix  A,  some  of  the  results 
given  in  [2]  need  to  be  modified.  The  improvement  factor  given 
by  (A. 15)  was  used  in  [2]  both  for  the  receiver  with  a  zero 
delay  suppression  filter  (whitening  filter  only)  and  for  the 
receiver  with  a  linear  phase  filter  (whitening  filter  followed 
by  its  matched  filter) ,  and  linear  phase  filter  based  on  the 
Welch  method.  Given  (A. 17),  however,  this  is  clearly  in¬ 
appropriate.  The  improvement  factors  for  these  latter  cases 
have  been  recalculated,  and  the  results  are  given  in  this 
appendix. 

Figure  B.l  is  the  newly  calculated  improvement  factor 
for  the  receiver  with  a  whitening  filter  followed  by  its 
matched  filter.  In  this  figure  the  horizontal  axis  is  the 
SNR/chip  without  filtering,  given  by 


*T 

T 


T 


This  figure  replaces  Figure  3.16  in  [2],  The  curves  in  Figure 
B.l  are  about  3dB  below  the  corresponding  curves  in  Figure 
3.16  of  [2].  This  implies  that  the  self-noise  term  in  the 
improvement  factor  is  in  this  case  the  dominant  noise  term, 
since  we  have  increased  the  self-noise  by  almost  3dB,  and  the 
resulting  change  in  the  improvement  factor  is  about  3dB. 


D  4  TRP  PREDICTION  FILTER 
0  15  TRP  PREDICTION  FILTER 

EXACT  fl. C. ,  100  TONES 


o?-0.01 


109 


Figure  3.41  in  [2]  will  also  be  affected  by  this  correction, 
although  we  have  not  replotted  it  here.  The  single  effect 
will  be  to  increase  the  two  self-noise  curves  in  this  figure 
by  3dB. 

Figure  B.2  is  a  figure  which  was  not  included  in  [2]. 
This  figure  contains  plots  of  the  improvement  ratio  for  both 
the  receiver  with  a  matched  filter  and  the  receiver  without 
the  matched  filter,  plotted  against  the  signal- to-white  noise 
ratio,  E{,/Nq. 

Figure  B.3  illustrates  the  improvement  factor  provided 
by  the  15-tap  linear  phase,  filter  obtained  using  the  Welch 
method.  This  figure  differs  little  from  the  corresponding 
figure  in  [2]  (Figure  3.25).  At  the  most,  there  is  1.5dB 
difference  in  the  curves.  This  is  because  the  self- noise 
term  is  not  the  dominant  term  with  this  filter  because  there 
is  significant  residual  interference  at  the  output. 
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